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ABSTRACT 

The  main  purpose  of  this  study  was  the  implementation  of  a  regional  flood  frequency 
analysis  for  Zion  National  Park  and  vicinity  areas.    Two  regional  models  were  implemented 
based  on  the  multivariate  lognormal  distribution  that  explicitly  accounts  for  the  space 
correlation  among  the  series  of  annual  peaks  in  the  region.    The  parameters  of  the  regional 
models  were  estimated  using  the  maximum  likelihood  procedure  based  on  the  observed  series 
of  peaks  from  a  subset  of  stations  that  are  known  to  preserve  the  assumed  structure  of  the 
regional  models.    The  models  have  the  capability  to  estimate  flood  quantiles  (maximum 
instantaneous  peaks)  at  ungaged  river  sites  based  on  readily  available  geographical  and 
meteorological  information.    Moreover,  the  models  can  also  be  used  to  increase  the 
reliability  of  estimated  quantiles  at  sites  with  very  short  records.    The  regional  models 
showed  very  good  fitting  capabilities  when  compared  with  empirical  frequency  distributions 
at  three  selected  stations  representative  of  catchments  at  low,  intermediate  and  high  elevation 
areas  of  the  Upper  Virgin  River  Basin. 
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1.0  INTRODUCTION 

1.1  Purpose  and  Scope  of  the  Study 

The  general  purpose  of  this  study  is  to  undertake  a  regional  peak  flows  frequency 
analysis  for  Zion  National  Park  and  its  surrounding  areas.    More  specifically,  the  objectives 
of  the  study  are  two  fold:  1)  to  estimate  peak  flow  quantiles  associated  to  given  exceedance 
probabilities  for  ungaged  river  sites,  and  2)  to  increase  the  reliability  of  estimated  peak  flow 
quantiles  for  sites  with  short  period  of  records. 

Zion  National  Park  is  located  in  South-western  Utah  where  surface  flow  records  are 
significantly  limited.    The  limitation  is  not  so  much  in  relation  to  the  number  of  gaging 
stations  (in  operation  and  discontinued),  but  regarding  the  length  of  the  peak  records 
available.    The  majority  of  the  gaging  stations  have  short  (less  than  50  years)  and  very  short 
records  (less  than  25  years).    Only  two  stations  in  the  Upper  Basin  of  the  Virgin  River  are 
classified  as  having  long  records  (more  than  50  years). 

Moreover,  the  reliability  of  estimated  quantiles  derived  from  single-site  analysis  is 
usually  low  because  of  the  limited  information  contained  in  a  single  series.    Customarily,  the 
problem  of  insufficient  information  at  a  specific  site  or  even  the  complete  lack  of  data  at  a 
site  has  been  tackled  by  using  additional  regional  information  from  neighboring  stations. 
Regional  flood  frequency  analysis  has  been  shown  to  be  very  effective  in  improving  the 
reliability  of  peak  flow  frequency  estimates  as  demonstrated  by  a  solid  body  of  publications 
in  the  field  that  starts  with  the  pioneering  works  by  Dalrymple  (1960)  and  Benson  (1962). 

The  methodology  for  flood  frequency  analysis  proposed  in  this  study  is  based  on  a 
multivariate  lognormal  model  that  explicitly  accounts  for  the  cross  correlation  among  the 
series  of  annual  peak  flows  in  the  region.    The  model  is  based  on  regional ization  assumptions 
derived  from  the  observed  series  of  peak  flows.    The  regional  structure  allows  the  model  to 
maintain  a  relatively  low  number  of  parameters  and  at  the  same  time  gain  in  reliability  by 
using  all  available  hydrological  data  in  the  region.    It  is  expected  that  the  limited  number  of 
years  of  data  for  most  streams  in  the  Zion  area  can  be  partially  compensated  by  this  regional 
approach. 


The  terms  'flood'  and  'peak  flow'  are  used  interchangeable  in  this  report.    However, 
we  would  like  to  make  a  distinction  between  peak  flows  with  intermediate  frequency  of 
occurrence  (once  every  2,  5,  10,  20  years)  versus  floods  of  devastating  effects  causing  severe 
human  and  economical  losses  (those  with  over  50  or  100  years  of  recurrence  interval).    This 
study  focuses  primarily  on  the  first  group  of  events.    This  is  in  accordance  to  the  general 
focus  of  other  geomorphological  and  ecological  studies  being  carried  out  in  the  Zion  area. 
Moreover,  there  is  a  general  agreement  among  experts  on  the  subject  that  suggest  that  the 
most  frequent  extreme  events  are  the  ones  that  have  the  highest  significance  on  the  stream 
ecosystem,  see  for  instance  Resh  et  al  (1988). 

This  report  includes  the  following  sections:    Section  1:    introduces  the  reader  to  the 
purpose  and  objectives  of  the  study  and  provides  a  brief  geographical  description  of  the  study 
area.    Section  2:    presents  the  basic  hydrological  information  used  in  the  study,  they  are 
series  of  maximum  instantaneous  peaks  and  regional  precipitation-frequency  data.    Section  3: 
provides  a  detailed  characterization  of  the  peaks  series  in  the  Zion  area.    It  also  explores 
some  criteria  for  hydrological  homogeneity  and  the  statistical  relationship  between  mean 
annual  peaks  and  hydro-meteorological  variables.    Section  4:    presents  the  formulation  and 
parameter  estimation  procedure  of  two  statistical  flood  frequency  models.    Section  5:    applies 
the  models  to  the  area  of  interest,  including  an  analysis  of  their  performance.    This  section 
also  provides  an  example  and  recommendations  for  using  the  models  to  estimate  flow  peak 
quantiles  at  ungaged  sites.    Finally,  Section  6:    gives  a  summary  of  the  study  and  final 
recommendations.    The  report  ends  with  literature  references  and  five  appendixes. 

1.2       Study  Area 

The  general  area  for  this  study  encompasses  a  large  portion  of  South-western  Utah, 
between  the  latitudes  N38°07'30"  and  N37°00'00\  and  the  longitudes  Wlllo15'00"  and 
W114°00'00",  covering  a  total  area  of  11,400  square  miles.    It  comprises  Washington 
County  almost  entirely,  and  parts  of  Kane,  Iron  and  Garfield  counties,  see  Figure  1.    The 
location  of  Zion  National  Park,  within  the  Upper  Virgin  River  Basin,  is  also  shown  in  Figure 
1.    Geographical  information  for  the  selected  watersheds  was  extracted  from  a  total  of  116 
topographic  maps  at  scale  1:24,000  (7.5-minutes)  prepared  by  the  U.  S.  Geological  Survey 
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(USGS).  An  inventory  of  the  topographic  maps  utilized  is  provided  in  Appendix  I.  Maps 
names  in  the  appendix  are  ordered  alphabetically,  including  also  their  corresponding  USGS 
code-reference  number. 

Our  study  region  encompasses  several  of  the  hydrologic  areas  of  the  southern  portion 
of  the  State  as  classified  by  the  Utah  Division  of  Water  Resources  (UDVVR,  1983,  pp.2). 
The  watersheds  of  interest  were  selected  based  on  the  availability  of  peak  flows  records  and 
some  criteria  of  hydrological  homogeneity  discussed  latter  in  Section  3.    The  following  lists 
the  state  hydrologic  areas  and  all  the  water  courses  contained  in  each  area  that  were  included 
in  the  study: 


Virgin  (9-7) 

-  East  Fork  of  the  Virgin  River 

-  North  Fork  of  the  Virgin  River 

-  North  Creek 

-  Virgin  River 

-  La  Verkin  Creek 

-  Ash  Creek 

-  South  Ash  Creek 


-  Crystal  Creek 

-  Deep  Creek 

-  East  Fork  Deep  Creek 

-  Mineral  Gulch 

-  Leeds  Creek 

-  Santa  Clara  River 

-  Moody  Wash 


Upper  Sevier  (5-6) 
-   Sevier  River 


Escalante  (8-4) 

-  Escalante  River 

-  Pine  Creek 

-  Deer  Creek 


-  Upper  Valley  Creek 

-  Twentymile  Wash 


Paria  (9-5) 

-  Paria  River 

-  Henrieville  Creek 

Kanab  (9-6) 

-  Kanab  Creek 

-  Johnson  Wash 


Cedar-Parowan  (6-2) 

-  Coal  Creek 

-  Center  Creek 

-  Summit  Creek 


Little  Creek 
Shirts  Creek 


>  Wahweep  (9-4) 

-   Coyote  Creek 

All  the  streams  listed  above  have  (active)  or  had  (discontinued)  USGS  surface  flow 
gaging  stations.    The  gaging  sites  provided  the  flows  records  used  in  the  regional  flood 
frequency  analysis.    Some  of  the  streams  have  flow  records  available  at  more  than  one 
location  as  indicated  later  in  Section  2.    A  list  of  the  USGS  topographic  maps  used  to 
delineate  and  study  each  of  the  subareas  depicted  in  Figure  1  is  included  at  the  end  of 
Appendix  I.    All  basic  geographical  and  hydrological  information  of  the  watersheds  and  the 
stream  network  is  provided  in  Section  2. 


2.0  PEAK  FLOWS  AND  PRECIPITATION  DATA 

2.1  Peaks  Flows  Data 

Streamflow  data  for  this  study  were  obtained  from  the  USGS  CD-Rom  database, 
named  Peak  Flows,  and  through  personal  communications  with  the  USGS  Water  Resources 
Division  in  Cedar  City,  Utah  (USGS,  1994).    Data  consist  of  maximum  annual  instantaneous 
flow  peaks,  expressed  in  cubic-feet  per  second  (cfs).    A  preliminary  selection  of  gaging  sites 
included  practically  all  records  available  within  the  study  area,  with  the  exception  of  just  a 
few  stations  having  flows  strongly  affected  by  human  activities.    Typically,  stations  located  at 
relatively  high  altitudes  have  short  records,  while  longer  records  are  likely  to  be  available  for 
the  larger  streams  having  large  drainage  areas.    This  is  also  the  case  in  the  Zion  region, 
where  periods  of  records  are  of  different  length  for  the  headwaters  than  for  the  lower  regions 
in  the  Upper  Basin  of  the  Virgin  River.    Nevertheless,  due  to  the  characteristics  of  the 
methodology  of  analysis  proposed  in  this  study,  more  is  to  be  gained  by  using  all  existing 
data  within  an  homogeneous  region  than  by  arbitrarily  discarding  stations  with  short  period 
of  records. 

2.1.1    Annual  Peak  Flows  Records 

Table  1  lists  39  gaging  stations  that  were  initially  considered  suitable  for  the  regional 
frequency  analysis.    The  first  column  in  Table  1,  site  no.,  provides  a  reference  number  to 
help  locate  the  gaging  stations  in  Figure  2,  which  shows  a  map  with  the  watersheds  of 
interest  and  the  main  drainage  network.    Station  type  indicates  whether  the  gaging  station 
was  operated  as  a  continuous-record,  denoted  by  'c\  or  only  as  a  crest-stage  partial  record 
station,  denoted  as  'p',  for  which  only  annual  peaks  are  published.    Latitude,  longitude  and 
gage  elevation  (datum)  information  correspond  to  the  site  where  the  gaging  station  is  located, 
obtained  from  USGS  records.    Drainage  area  is  listed  next.    Some  minor  differences  in 
watersheds  areas  were  detected  between  the  values  provided  by  the  USGS  and  those 
computed  by  the  Geographic  Information  System  (GIS)  developed  during  this  study.    Base 
flow  levels,  established  by  the  USGS  to  create  partial-duration  series  of  peak  flows,  are 
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provided  in  Table  1  only  for  those  stations  with  long  and  continuous  records.    Period  of 
record  provides  the  first  and  last  water-year  of  available  data,  with    Y1  and  'N'  indicating 
whether  a  series  is  continuous  or  discontinuous  during  the  period,  respectively.    The  number 
under  Count  indicates  the  length  of  the  peak  annual  series  in  years.    Notice  that  only  15  out 
of  the  39  stations  are  still  reporting  data  during  water-year  1993.    For  easy  visualization,  a 
bar  graph  showing  graphically  the  length  of  usable  records  is  presented  in  Figure  3. 

Table  1  also  designates  the  flood  region  where  each  catchment  with  a  gaged  site  is  in. 
This  is  based  on  the  Thomas  and  Lindskov  (1983)  classification  of  flood  regions  for  the  State 
of  Utah.    To  determine  the  flood  region,  the  user  needs  to  estimate  the  datum  and  mean 
basin  elevation  of  the  site  of  interest.    The  datum  is  the  elevation  at  the  flow  gaging  site,  and 
the  mean  basin  elevation  is  determined  for  the  drainage  basin  which  contributes  streamflow 
to  the  site.    Approximate  boundaries  of  the  flood  regions  intersecting  the  study  area  are 
shown  in  Figure  2.    Further  discussion  of  this  subject  is  provided  in  Section  3.1.    Finally,  a 
few  remarks  are  included  in  Table  1  concerning  the  degree  of  water  diversion  and/or 
regulation  in  the  streams.    The  Upper  Basin  of  the  Virgin  River  and  adjacent  watersheds 
have  not  been  heavily  impacted  by  human  activities;  particularly  in  relation  to  flood  control 
structures  which  can  drastically  reduce  the  magnitude  of  flow  peaks.    However,  the  lack  of 
accurate  and  systematic  water  diversion  records  prevented  us  from  analyzing  the  effect  of 
water  diversions  in  the  magnitude  of  flow  peaks.    Unfortunately,  the  months  with  the  largest 
demand  for  irrigation  water,  June  through  September,  coincide  with  the  period  when  flow 
peaks  reach  their  maximum  frequency  and  intensity,  July  through  September.    All  available 
peak  flows  data  were  tabulated  and  included  in  the  report  as  Appendix  II  (see  values  in 
normal  type  face).    The  first  group  of  basic  statistics  reported  at  the  bottom  of  the  page 
corresponds  to  the  series  of  measured  peaks. 

2.1.2    Extension  of  Annual  Peak  Flows  Records 

The  validation  of  some  of  the  underlying  assumptions  of  the  multivariate  model  to  be 
applied  here  requires  a  preliminary  estimation  of  the  long-term  mean  annual  peak  (MAP)  at 
each  gaging  site.    The  difficulty  in  estimating  MAPs  arises  from  the  unequal  length  of  the 
records  available.    As  mentioned  before,  few  records  extend  over  60  or  70  years,  but  the 
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great  majority  are  considered  short  or  very  short  records.    Prior  to  the  estimation  of  MAPs, 
it  is  required  that  some  records  be  extended  to  cover  a  common  base-period.    This  approach 
makes  data  more  comparable  by  diminishing  the  effect  that  natural  cycles  in  weather 
conditions  may  have  in  those  stations  with  short  records.    This  is  a  standard  procedure 
utilized  by  the  USGS  when  applying  index-flood  methods  (Dalrymple,  1960).    The  base-  or 
common-time  period  should  be  chosen  as  long  as  possible.    The  data  set  available  for  the 
region  at  hand  dictates  the  maximum  possible  length  of  the  base-period.    For  instance, 
observing  the  bar  graph  in  Figure  3  it  is  concluded  that  the  base  period  should  extend  from 
1959  to  1993,  since  the  largest  number  of  records  start  in  1959. 

Simple  and  multiple  linear  regression  techniques  allowed  us  to  fill  in  missing  peaks  at 
several  stations.    Figure  4  shows  a  linear  correlation  between  annual  peaks  at  two  gaging 
sites  along  the  Paria  River.    The  correlative  association  between  the  two  series  is  typical  for 
the  study  area  when  gaging  sites  are  relatively  close  or  along  the  same  water  course.    For 
pairs  of  stations  separated  a  large  distance  the  linear  correlation  of  peaks  deteriorates  rapidly. 
In  some  cases  linear  correlations  could  not  be  detected.    This  is  due  to  the  random 
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occurrence  in  time  and  space  of  storms  producing  flow  peaks,  and  the  low  probability  of  an 
intense  summer  storm  hitting  a  particular  site,  particularly  in  small  catchments.    Given  the 
relative  low  values  of  correlation  coefficients  between  annual  peak  series,  visual  inspection 
became  an  important  criterion  to  help  determine  whether  a  correlative  association  between 
two  stations  was  acceptable  or  not  for  the  purpose  of  record  extension. 

The  correlation  models  finally  adopted  included  random  components  to  help 
reproducing  the  full  variance  of  the  dependent  variable.    With  only  a  few  exceptions,  actual 
peak  measurements  were  used  to  fill  in  data  at  another  station,  avoiding  the  use  of  previously 
estimated  peaks  for  further  record  extension.    Whenever  possible,  records  were  also  extended 
or  filled  in  prior  to  1959.    This  was  done  only  for  stations  that  showed  significant  linear 
correlation  with  other  existing  records  during  the  period  1910-1959.    The  extended  data  base 
of  peak  flows  is  also  presented  in  Appendix  II.    Filled  in  values  are  highlighted.    The 
original  data  base  was  extended  from  a  total  of  847  year-station  values  to  1249  values,  a 
47.4%  increase.    Basic  statistics  for  the  extended  series  are  also  included  in  Appendix  II  at 
the  bottom  of  the  page,  below  the  set  of  statistics  for  the  systematic  series. 

2.2       Precipitation-Frequency  Data 

Since  precipitation  is  the  primary  cause  of  peak  flows,  it  is  imperative  that  some  form 
of  precipitation  be  part  of  the  regional  flood  frequency  analysis.    In  New  England,  where 
rainfall  is  fairly  uniform,  Benson  (1962)  found  that  rainfall  intensity,  though  a  statistically 
significant  variable,  was  not  one  of  the  most  important  factors  affecting  flow  peaks.    On  the 
other  hand,  in  the  southwest  where  as  Benson  (1964)  says,  "...  the  outstanding  characteristic 
of  precipitation  is  its  variability",  precipitation  is  expected  to  be  an  important  parameter  in 
the  determination  of  mean  annual  flow  peaks  in  the  region. 

Heavy  summer  rains  in  southern  Utah  create  the  peak  flows  with  the  largest  intensity 
in  the  region.    According  to  Woolley  (1946),  there  are  two  distinct  types  of  storms:  (1)  the 
steady  or  cyclonic  storms,  which  are  not  very  frequent  and  may  result  in  a  large  total 
precipitation  but  are  spread  over  a  protracted  period  and  have  relatively  small  intensities;  and 
(2)  thunderstorms,  which  generally  include  the  violent  downpours  of  brief  duration  (20  to  30 
minutes)  and  high  intensity  commonly  known  as  cloudbursts.    The  great  majority  of  the 
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cloudburst  events  occur  in  July  and  August,  the  warmest  months,  while  a  few  events  occur  in 
June  and  September.    Flow  peaks  induced  by  cloudbursts  can  exceed  in  magnitude  the 
normal  spring  freshets.    They  are  characterized  by  very  flashy  peaks  followed  by  quickly 
diminishing  flows.    Most  of  the  cloudburst  peaks  in  the  Zion  region  originate  below  the 
8,000  ft  altitude,  and  mainly  in  streams  exposed  to  the  direction  of  storm  travel.    The  very 
local  effect  of  high  intensity  summer  thunderstorms  can  yield  erroneous  results  when  flood 
frequency  analysis  is  carried  out  at  a  single  site,  particularly  if  dealing  with  small  drainages 
and  short  periods  of  observation. 

There  are  several  indexes  of  precipitation  that  can  be  used  in  the  study  regarding  the 
distribution  of  mean  annual  peaks  in  a  region.    Depending  on  the  meteorological  mechanism 
generating  the  peaks,  different  precipitation  indexes  can  be  of  use.    For  instance,  for  snow 
melt-flood  areas,  mean  annual  precipitation  has  been  found  significantly  correlated  with 
annual  peaks  (Benson,  1964).    On  the  other  hand,  for  rainfall-induced  peaks,  a  precipitation 
index  involving  some  form  of  frequency  and  duration  should  be  highly  correlated  with  the 
instantaneous  peak  discharge  of  a  given  frequency.    This  study  considers  two  indexes  from 
the  latter  class,  the  2yr-24hr  precipitation  and  the  50%-24hr  precipitation.    Maps  of 
isopluvial  lines  for  both  indexes  were  obtained  from  the  National  Weather  Service  (NOAA, 
1973a).    Precipitation  maps  presently  available  for  the  State  of  Utah  were  generated  by  the 
NWS  from  records  compiled  from  1910  to  1969.    At  present,  the  NWS  is  in  the  process  of 
updating  the  precipitation-frequency  maps  including  information  from  the  last  two  decades. 
Isopluvial  lines  were  digitized  and  incorporated  as  an  additional  thematic  layer  into  the 
regional  GIS.    A  reproduction  of  one  of  the  precipitation  layers  is  shown  in  Figure  5. 

The  two  adopted  indexes  represent  different  precipitation  conditions.    The  2yr-24hr 
precipitation  map  was  generated  based  on  precipitation  data  from  all  months  of  the  year, 
without  regard  to  the  type  of  precipitation:  rain,  rain  and  snow  mixed  or  all  snow.    The  2 
year  return  period  for  precipitation  was  chosen  to  match  the  expected  recurrence  interval  of 
mean  annual  peaks,  which  in  general  is  estimated  to  lie  between  2  to  3  years.    On  the  other 
hand,  the  0.50  probability-6hr  precipitation  map  considers  only  precipitation  events  from  the 
months  of  May  to  October.    This  is  equivalent  to  using  a  data  series  made  up  of  precipitation 
events  that  are  exclusively  rain  (rainfall  only  series).    Because  summer  thunderstorms  have 
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relatively  short  duration,  the  precipitation  map  with  the  shortest  duration  available  for  this 
series  (6-hr)  was  chosen. 

The  importance  of  snow  in  precipitation-frequency  for  the  region  will  dictate  which  of 
the  two  indexes  will  apply  better  for  the  purpose  of  the  study.    It  is  known  that  at  very  high 
elevations  snow  can  account  for  over  50%  of  the  normal  annual  precipitation.  However,  that 
doesn't  imply  that  snow  contributes  in  the  same  proportion  to  the  annual  series  of  maximum 
6hr  or  24hr  precipitation  amounts.    A  study  carried  out  by  the  NWS  (NOAA,  1973b)  for  the 
State  of  Utah  and  other  Rocky  Mountain  States  analyzed  the  effect  of  snow  occurrences  on 
the  mean  annual  precipitation.    Precipitation-frequency  analysis  derived  from  two  data  sets 
were  compared.    One  series  containing  precipitation  amounts  for  each  year,  regardless  of  the 
type  of  precipitation,  and  a  second  series  restricted  solely  to  rainfall  events.    The  study 
concluded  that,  for  2yr  return  period  and  durations  of  24hr  or  less,  the  contribution  of  snow 
amount  to  precipitation-frequency  values  in  the  rainfall  only  series  averages  10  to  15% 
smaller  than  the  series  using  all  precipitation  amounts.    Same  analysis  for  the  25yr  return 
period  produced  differences  somewhat  smaller.    This  indicates  a  tendency  for  the  rare  short- 
duration  storm  to  occur  in  the  warmer  months  as  rain. 
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3.0  REGIONAL  CHARACTERISTICS  OF  PEAK  FLOWS  SERIES 

3.1  Hydrometeorological  Mechanisms  Producing  Peak  Flows 

A  critical  aspect  of  flood  frequency  analysis  is  the  understanding  of  the 
hydrometeorological  mechanisms  producing  the  peaks.    Numerous  studies  in  the  past  were 
focused  in  that  direction.    Specifically  for  the  State  of  Utah,  Thomas  and  Lindskov  (1983) 
pointed  out  that  annual  peak  flows  at  many  streams  in  Utah  can  result  from  rainfall,  snow 
melt  or  rain  on  snow,  and  consequently,  a  statistical  flood  frequency  relation  determined 
from  an  array  of  annual  peaks  from  a  mixed  population  usually  does  not  adequately  fit  the 
larger  peak  flows.    The  two  authors  carried  out  an  extensive  study  with  the  purpose  of 
defining  flood  regions  with  unique  statistical  characteristics.    They  concluded  that  the  upper 
altitudinal  limit  for  floods  caused  by  intense  thunderstorms  increased  from  North  to  South  in 
the  State.    Flood  regions  were  defined  based  on  residual  patterns,  mean  basin  elevation  of  the 
gaging  stations,  and  on  the  type  of  floods  that  have  occurred  at  the  gaging  sites  (snow  melt, 
thunderstorm,  frontal  rainfall,  or  combination  thereof).     The  reader  can  find  a  complete  map 
of  the  State  of  Utah  showing  all  flood  regions  in  page  12  of  the  above  reference.    In  an 
earlier  work,  Wooley  (1946,  p. 50)  also  indicated  that  "..  Most  of  the  cloudburst  floods 
(produced  by  high  intensity,  short  duration  storms)  originate  in  regions  below  the  8,000  ft 
altitude,  in  streams  more  or  less  directly  exposed  to  the  direction  of  storm  travel". 

Boundaries  of  flood  regions  intersecting  our  study  area  were  approximately 
reproduced  in  Figure  2.    Three  flood  regions  cross  the  area,  they  are:  the  Low  Plateau 
region  (LP),  the  High  Plateau  region  (HP)  and  the  High  Elevation  subregion.    Based  on 
Thomas  and  Lindskov's  classification  criteria,  only  4  gaging  sites  are  part  of  the  HP  region, 
whereas  the  remaining  35  sites  are  contained  within  the  LP  region,  see  Table  1.    The  four 
HP  sites  are  located  in  the  North  Fork  of  the  Virgin  River  subbasin,  upstream  from  the  Park, 
next  to  the  watershed  divide  in  the  Cascades  Spring  area.    Catchments  within  the  HP  region, 
south  of  38°  latitude,  should  have  a  mean  basin  elevation  greater  than  8,000  ft  and  a  study 
site  datum  greater  than  7,000  ft.    This  is  a  mixed  population  flood  area  where  infrequent 
thunderstorm  floods  can  have  greater  magnitude  than  snow  melt  floods.    Catchments  within 
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the  LP  region,  also  south  of  38°  latitude,  must  have  a  mean  basin  elevation  less  than  8,000  ft 
and  a  study  site  datum  less  than  7,000  feet.    Summer  thunderstorms  produce  most  of  the 
large  magnitude  floods.    Snow  melt  floods  are  rare  and  usually  small. 

Thomas  and  Lindskov  (1983)  flood  areas  can  be  interpreted  as  a  stratified 
classification  of  the  region  where  different  flood  producing  mechanism  are  dominant.    For 
instance,  for  basins  south  of  latitude  39°,  flood  frequency  relation  at  elevations  below  7,000 
feet  are  dominated  by  thunderstorms.    Above  9,000  feet,  flood  frequency  relations  will  be 
controlled  by  snow  melt,  except  for  recurrence  intervals  larger  than  100  years  where 
thunderstorms  become  significant  again.    The  intermediate  zone,  between  7,000  feet  and 
9,000  feet,  will  have  a  mixed  population.    Another  interesting  conclusion  from  Thomas 
(1985),  who  proposes  a  composite  flood-frequency  analysis  for  gaging  sites  with  mixed- 
populations  of  peak  flows,  is  that  for  basins  with  gages  around  6,000  ft  in  elevation  and 
below  (similar  to  sites  of  interest  in  Zion  National  Park),  flood-frequency  relations  derived 
from  the  annual  maximum  peaks  (without  a  prior  disaggregation  of  populations)  is  adequate 
for  the  commonly  used  range  of  recurrence  interval  of  10  to  100  yrs. 

Jarrett  (1993)  also  found  latitude-dependent  elevation  limits  for  different 
meteorological  mechanisms  producing  floods  in  the  Rocky  Mountains  region.    Pons  (1992) 
was  also  concerned  with  peak  flows  populations  of  mixed-origin.    During  his  investigation, 
Pons  selected  four  areas  of  the  country  with  uniform  hydrometeorological  conditions  to  avoid 
problems  associated  with  multiple  flood  producing  mechanisms.    The  segregation  of  flow 
records  in  separate  populations  to  reflect  single  peak  producing  mechanisms  is  certainly  an 
ideal  and  desirable  condition.    However,  the  scarcity  of  peak  flows  records  available  for 
South-Western  Utah  precludes  that  approach  for  the  majority  of  the  gaging  sites. 

3.2       Seasonal  Distribution  and  Duration  of  Peak  Flows 

Although  Thomas  and  Lindskov  (1983)  when  carrying  out  their  multiregression 
analysis  did  not  explicitly  differentiated  between  snow  melt  and  rainflood  peaks   (as  done  for 
instance  by  Peterson  (1985)  who  estimates  snow  melt  and  thunderstorm  flood  magnitudes  for 
specific  drainages),  they  have  implicitly  taken  into  account  the  nature  of  the  peak  flows  by 
incorporating  mean  basin  elevation  (together  with  drainage  area)  in  the  sets  of  regional 
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regression  equations.    The  regional  classification  procedure  implemented  by  Thomas  and 
Lindskov  was  adopted  in  this  study  as  a  preliminary  criterion  for  selecting  gaging  sites.    It 
should  provide  some  degree  of  hydrological  homogeneity  among  peak  flows  records  within 
the  study  area. 

The  objective  of  this  section  is  to  investigate  the  hypothesis  that  peaks  within  the 
study  area  are  derived  from  a  unique  population.    Although  this  is  not  a  basic  requirement  of 
the  proposed  methodology  of  analysis,  it  helps  remove  sources  of  uncertainty  concerning  the 
proper  fitting  of  the  regional  probabilistic  model.    Two  representative  watersheds  with  long 
and  well  established  records  were  selected  for  comparison.    The  first  watershed  corresponds 
to  the  North  Fork  of  the  Virgin  River.    This  catchment  has  a  large  drainage  area  extending 
from  high  elevation  headwaters  to  a  low  elevation  gaging  site  at  3,970  ft  near  the  town  of 
Springdale  (site  23,  ID:  09405500).    The  drainage  network  is  mainly  oriented  in  the  north- 
south  direction  and  subject  to  a  strong  gradient  of  total  precipitation  that  ranges  from  29 
inches  at  the  headwaters  to  15  inches  at  the  outlet  per  year  (see  Figure  11  in  Diaz  and 
Ismael,  1993).    By  contrast,  the  second  watershed  was  chosen  among  the  ones  with  the 
highest  elevation,  and  small  drainage  area.    The  selected  watershed  is  Coal  Creek  near  Cedar 
City  (site  38,  ID:  10242000).    It  has  an  area  of  81  mi2,  with  its  outlet  at  6,000  ft,  and  a 
mean  basin  elevation  of  8,640  ft.    The  location  of  the  watersheds  can  be  found  in  Figure  2. 

Figure  6  displays  the  relative  magnitude  and  the  time  of  occurrence  of  flow  peaks 
registered  at  sites  23  and  38,  top  and  bottom  graphs  respectively.    The  horizontal  axis 
indicates  the  month  of  the  Julian-year  at  which  a  flow  peak  occurred,  showing  the  empirical 
distribution  of  seasonal  peak  occurrence.    The  vertical  axis  plots  the  ratio  between  the 
magnitude  of  an  instantaneous  peak  to  the  corresponding  mean-daily  flow.    This  graph 
provides  information  about  the  degree  of  "flashiness"  of  the  peak  events  and  their  time  of 
occurrence.    There  is  a  clear  tendency  for  the  peaks  to  occur  during  the  summer  season,  July 
through  September  at  both  watersheds,  though  peaks  also  occur  year  around  mostly  at  low 
elevation  basins  (see  top  graph).    Furthermore,  visual  inspection  of  the  annual  hydrographs  at 
Springdale  allowed  us  to  ascertain  the  preponderance  of  peak  flows  of  rainfall  origin  over 
those  of  snow  melt  origin.    It  should  be  mentioned  that  partial  series  of  peaks,  instead  of 
annual  series,  were  used  for  creating  Figure  6.    Partial  series  include  annual  series  plus 
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Fig. 6   Relative  Magnitude  and  Time  of  Occurrence  of  Peak  Flows 
(top)   N.  Fk.  of  the  Virgin  River  near  Springdale,  and 
(bottom)   Coal  Creek  near  Cedar 
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annual  secondary  peaks  above  a  certain  threshold  flow  (Qb=  550  cfs  at  site  38  and  Qb«  820 
cfs  at  site  23).    The  base  flow  level  is  usually  selected  such  as  on  the  average  about  three 
secondary  peaks  per  year  are  included.    When  comparing  the  results  of  the  analysis  from 
both  types  of  series  it  was  found  that  practically  the  same  magnitude  of  ratios  and 
frequencies  were  depicted.    Hence,  partial  series  were  finally  adopted  because  of  the  larger 
number  of  peaks  available,  yielding  more  reliable  estimates  of  the  relative  magnitude  and 
time  occurrence  of  peaks. 

The  peak-daily  flow  ratio  can  be  interpreted  as  an  indicator  of  the  duration  of  the 
peaks  (flashiness).    The  larger  the  ratio  the  shorter  the  duration  of  the  peak,  and  vice  versa. 
During  the  summer  months,  the  ratio  reaches  magnitudes  much  larger  than  during  the  rest  of 
the  year.    High  intensity  and  short  duration  summer  storms  give  rise  to  a  sharp  increase  and 
fast  decay  of  the  hydrographs,  with  peaks  reaching  up  to  about  35  times  the  mean-daily  flow 
for  the  same  day  at  Springdale,  and  up  to  about  60  times  at  Coal  Creek.    On  the  other  hand, 
during  the  snow  melt  season,  usually  occurring  from  late  March  through  June,  instantaneous 
peaks  rarely  exceed  mean-daily  flows  more  than  4  times.    Remember  that  mean-daily  flows 
during  the  snow  melt  season  are  substantially  larger  than  during  the  rest  of  the  year.    This 
factor  is  mostly  responsible  for  the  low  values  of  the  ratios.    In  turn,  peak  events  with  the 
largest  ratios  within  the  snowfall  or  snow  melt  seasons  correspond  mostly  to  rainfall  on  snow 
events,  whereas  peak  events  of  exclusively  snow  melt  origin  are  the  fewest  and  have  the 
smallest  ratios  (close  to  1).    The  repetitive  nature  of  the  snow  melt  peaks  is  confirmed  by  the 
low  variance  of  the  ratios  during  the  months  of  April  to  June.    Conversely,  summer  peaks 
display  the  largest  variance  as  a  result  of  the  random  nature  of  the  rainfall  events. 

The  peak-flow  ratio  can  also  be  used  to  obtain  indirect  estimates  of  the  expected 
magnitude  of  peaks  for  each  season.    Consider  for  instance  the  month  of  August  at 
Springdale.    The  peak-flow  ratio  for  that  station  during  the  8th  month  ranges  from  1.7  to 
36.2,  with  an  average  value  of  9.6.    From  a  previous  study  we  know  that  base  flow  during 
August  has  a  relatively  uniform  value  of  64  cfs.    By  multiplying  the  last  two  figures,  a  peak 
flow  of  614  cfs  is  obtained.    This  constitutes  a  reasonable  estimate  for  an  average  rainfall 
peak  to  be  expected  during  the  month  of  August.    Similar  computation  for  the  month  of  May 
yields  an  average  peak  flow  of  500  cfs. 
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In  summary,  peak  flows  in  the  Zion  region  display  mixed-population  characteristics 
resulting  from  snow  melt  and  rainfall  events,  with  the  latter  as  the  predominant  mechanism. 
The  seasonal  distribution  and  duration  of  peak  flows  at  two  stations  representative  of  high 
and  intermediate  elevations  also  show  similar  characteristics.    Peak-daily  flow  ratios  are  low 
and  very  constant  (small  variance)  during  the  snow  melt  season,  and  large  and  highly 
variable  during  the  months  with  convective  summer  storms  (cloudbursts).    There  is  little 
evidence  of  the  occurrence  of  peaks  outside  the  period  May-September  in  high  elevation 
basins  (with  datum  over  6,000  ft).    However,  as  mean  basin  elevation  decreases,  in  addition 
to  snow  melt  peaks  and  summer  thunderstorms  peaks,  flow  peaks  are  also  induced  by  frontal 
systems  moving  across  the  state  during  the  snow  accumulation  period. 

3.3       Seasonal  Distribution  of  the  Number  of  Peaks 

Specific  environmental  studies  require  information  about  the  probability  distribution  of 
the  number  of  peaks  occurring  at  different  periods  of  the  year.    For  this  purpose,  the  effect 
of  seasonality  on  the  occurrence  of  flow  peaks  and  the  underlying  probability  law  have  been 
examined.    Zelenhasic  (1970)  presented  a  theoretical  approach  for  the  analysis  of  the 
stochastic  properties  of  flow  peaks.    This  approach  considers  the  number  of  peaks  occurring 
in  a  given  time  interval  as  well  as  the  time  when  the  peaks  occur  as  random  variables. 
Furthermore,  it  demonstrates  that  the  probability  distribution  governing  the  occurrence  of 
peaks  follows  a  Poisson  process. 

The  above  cited  methodology  has  been  applied  in  this  study  to  gain  further  insight 
into  the  probabilistic  structure  of  the  seasonal  occurrence  of  peak  flows  in  the  Upper  Virgin 
River  basin.    The  gaging  station  chosen  for  analysis  is  Virgin  River  at  Virgin  (Site  25,  ID: 
09406000).    This  station  measures  flow  peaks  from  the  entire  Upper  Virgin  River  basin  (956 
mi2).    It  has  the  longest  flow  record  available,  including  a  relatively  large  number  of 
secondary  peaks  in  record.    The  data  are  provided  in  Appendix  III.    The  Virgin  River  data 
shows  an  adequate  record  of  peaks  flows,  except  during  the  early  period  1910  to  1951,  when 
secondary  peaks  were  not  tabulated.    For  the  computations,  the  water-year  is  divided  into 
nine  periods,  eight  periods  40  days  long  and  one  period  of  45  days.    The  beginning  and  end 
of  the  periods  were  selected  to  accommodate  the  observed  seasonal  distribution  of  peaks  in 
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the  Virgin  River  basin.    Observed  frequencies  of  peaks  are  obtained  from  the  partial  series 
for  all  nine  periods.    Table  2  provides  a  summary  of  the  computations. 

For  example,  water-year  1961  contains  the  third  largest  annual  peak  on  record 
(13,500  cfs),  and  the  largest  number  of  secondary  peaks  registered  (9  peaks).    Six  out  this 
nine  secondary  peaks  correspond  to  the  same  period.    This  is  reflected  in  Table  2  (see  Period 
7)  as  one  case  observed  (Fob=l)  when  six  peaks  (k  =  6)  occurred  during  the  same  period.    In 
addition  to  the  observed  frequency,  Table  2  provides  the  theoretical  Poisson  distribution  (FJ 
of  peak  occurrences  during  each  period.    The  Poisson  distribution  has  a  single  parameter  k 

Table  2.    Observed  and  Fitted  Poisson  Distribution  of  the  Number  of  Peaks  for  Nine  Periods 


Per:  1 

Dec  4 -Jan  17 

k 

Fob 

Fth 

0 

144 

144  083 

1 

5 

4.835 

2 

0 

0.081 

3 

0 

0.001 

4 

0 

0000 

5 

0 

6 

0 

Sum: 

149 

149.000 

mean: 

0.034 

variance: 

0.033 

lambda: 

000075 

Per:  2 

Jan  18- 

reb26 

k 

Fob 

Fth 

0 

145 

143  158 

1 

2 

5.726 

2 

2 

0.115 

3 

0 

0.002 

4 

0 

0.000 

5 

0 

6 

0 

Sum: 

149 

149.000 

mean: 

0.040 

variance: 

0.066 

lambda: 

0.00101 

Per:  3 

Feb  27  - 

Apr  7 

k 

Fob 

Fth 

0 

139 

138.372 

1 

9 

10.240 

2 

1 

0.379 

3 

0 

0.009 

4 

0 

0000 

5 

0 

6 

0 

Sum: 

149 

149.000 

mean: 

0.074 

variance: 

0.082 

lambda: 

0.00185 

Per:  4 

Apr  8 -May  17 

k 

Fob 

Fth 

0 

137 

132.946 

1 

8 

15.156 

2 

3 

0.864 

3 

1 

0.033 

4 

0 

0.001 

5 

0 

0000 

6 

0 

Sum: 

149 

149.000 

mean: 

0.114 

variance: 

0.183 

lambda 

0.00285 

Per:  5 

May  18- 

Jun25 

k 

Fob 

Fth 

0 

144 

144.083 

1 

5 

4.835 

2 

0 

0.081 

3 

0 

0.001 

4 

0 

0.000 

5 

0 

6 

0 

Sum: 

149 

149.000 

mean: 

0.034 

variance: 

0.033 

lambda: 

0.00084 

Per:  6 

Jun  26  - 

<Vug  5 

k 

Fob 

Fth 

0 

123 

120.203 

1 

21 

25.815 

2 

4 

2.772 

3 

1 

0.198 

4 

0 

0.011 

5 

0 

0.000 

6 

0 

Sum: 

149 

149.000 

mean: 

0215 

variance: 

0.264 

lambda: 

0.00537 

Per:  7 

Aug  6  -  Sep  14 

k 

Fob 

Fth 

0 

114 

109.423 

1 

30 

33.782 

2 

2 

5.215 

3 

2 

0.537 

4 

0 

0.041 

5 

0 

0.000 

6 

1 

Sum: 

149 

148.997 

mean: 

0.309 

variance: 

0.526 

lambda: 

0.00772 

Per:  8 

Sep  15- 

Oct  24 

k 

Fob 

Fth 

0 

136 

136.551 

1 

13 

11.914 

2 

0 

0.520 

3 

0 

0.015 

4 

0 

0.000 

5 

0 

6 

0 

Sum: 

149 

149.000 

mean: 

0.087 

variance: 

0.080 

lambda: 

0.00218 

Per:  9 

Oct  25  -  Dec  3 

k 

Fob 

Fth 

0 

139 

138  396 

1 

9 

10.217 

2 

1 

0.377 

3 

0 

0009 

4 

0 

0  000 

5 

0 

6 

0 

Sum: 

149 

149.000 

mean: 

0.074 

variance: 

0.082 

lambda: 

0.00185 

22 


(lambda).    Its  mean  and  variance  are  both  equal  to  this  parameter.    The  parameter  X  is 
estimated  for  each  period  by  the  method  of  moments  using  information  from  the  available 
samples.    The  conformance  of  the  theoretical  to  the  observed  distributions  was  verified  by 
the  x2-test.    The  computed  x2-values  were  below  the  critical  x2-values  at  the  5  %  level  of 
significance  for  each  of  the  nine  periods  considered.    Sample  means  and  sample  variances  for 
the  number  of  occurred  peaks  during  the  nine  periods  of  the  year  are  given  in  Table  2.    As 
expected,  similar  values  for  the  two  parameters  are  obtained  at  each  period.    The  nine 
parameters  k ,  expressed  as  number  of  events  per  day,  are  plotted  in  Figure  7  as  the 
abscissas  at  the  center  of  each  time  period.    The  A,  values  are  then  fitted  by  Fourier  Series  to 
generate  a  continuous  function.    The  fitting  yields  a  deterministic  periodic  function  of  time, 
X  ( t)  ,  with  a  one-year  period  as  shown  in  Figure  7. 
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Zelenhasic  (1970)  expresses  the  expected  value  of  the  number  of  exceedances  r\  ( t] 
during  a  time  interval  [0,t]  as, 
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B[i?(t)]    = 


X(s)    ds  -  A(t)  (3-1) 


where  A(  t)  depends  on  A.  (s)  which  is  the  periodic  nonlinear  function  portrayed  in  Figure  7. 
Then,  the  probability  that  exactly  k  peaks  occur  in  a  fixed  period  [0,t]  is  estimated  by, 

P[K=k]    =     [A(t)]*exp[-A(t)]  f  k,012  M  (3.2) 

Equation  (3.2)  can  also  be  used  to  estimate  the  probability  that  k  peaks  will  occur  in 
a  fixed  time  interval  [t^tj,  where  t,  and  t2  correspond  for  instance  to  the  first  and  last  day  of 
a  given  month.    As  an  illustration,  we  will  determine  the  probability  that  at  least  two  peaks 
(two  or  more  peaks)  will  occur  during  the  month  of  August  ,  the  month  with  the  largest 
occurrence  of  peaks.    From  Eq.(3.1)  and  for  the  time  interval  corresponding  to  the  month  of 
August  (indicated  in  Figure  7  with  a  dotted  pattern),  the  average  rate  (of  arrival)  of  the 
Poisson  process  is  AAug=  0.25.    Then, 


P[K>2]    =  1-P[K<1]    =  1-FK{1, 


k=l 


=   1  -     k=0 


Y,  e-°-250.25* 


(3.3) 


k\ 

=  0.027 

indicates  that  there  is  only  a  2.7%  of  probability  that  two  or  more  peaks  above  the  base 
discharge  (1,600  cfs)  will  occur  at  Virgin  during  August.    Similarly,  the  probability  of 
having  at  least  one  peak  is  equal  to  22%,  and  the  probability  of  not  having  any  peak  during 
August  is  78%.    For  the  month  with  the  lowest  number  of  peaks,  January,  similar 
computations  yield:  AJan=  0.02,  a  negligible  probability  for  the  case  of  two  or  more  peaks, 
and  2%  and  98%  for  one  and  zero  peaks,  respectively. 
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3.4       Graphical  Estimation  of  Mean  Annual  Peaks 

As  mentioned  briefly  in  Section  2.1.2,  estimates  of  the  long-term  mean  annual  peak 
(MAP)  at  each  gaging  site  should  be  determined  in  order  to  assess  the  conditions  of 
hydrological  homogeneity  in  the  study  area.    This  task  is  performed  prior  to  the  application 
of  the  multivariate  flood  frequency  model.    Peaks  records  were  first  extended  to  cover  the 
longest  possible  common  base-period.    Then,  each  extended  series  of  annual  peaks  were 
subject  to  a  preliminary  single-site  frequency  analysis.    The  preliminary  analysis  consisted  on 
plotting  annual  peaks  discharges  versus  recurrence  intervals  for  each  gaging  station. 
Recurrence  intervals  were  computed  by  arranging  the  data  in  decreasing  order  of  magnitude 
and  assigning  a  plotting  position  to  the  ranked  values.    The  standard  Weibull  formula  was 
used   for  this  purpose.    Although  estimated  peaks  (obtained  using  the  regression  procedure) 
are  part  of  the  annual  series  used  to  rank  the  flows,  they  are  not  included  in  the  frequency 
plots.    Plots  of  the  single-site  frequency  curves  are  included  in  this  report  as  Appendix  IV. 

Using  the  plots  in  Appendix  IV,  the  MAP  at  each  gaging  site  is  estimated  graphically 
by  determining  the  intersection  of  a  visually  fitted  frequency  curve  with  the  vertical  line 
corresponding  to  a  recurrence  interval  of  2.7  yr.    This  recurrence  interval  approximately  to 
the  recurrence  of  the  mean  for  the  lognormal  and  log-Pearson  type-Ill  distributions. 
Frequency  points  were  carefully  fitted  in  the  1  yr  to  10  yr  recurrence  interval.    This 
graphical  procedure  to  estimate  the  mean  annual  flood  gives  more  weight  to  peaks  occurring 
in  the  neighborhood  of  the  mean.    This  approach  contrasts  with  a  simple  arithmetic  average 
of  the  peaks  which  would  be  heavily  biased  by  infrequent  peaks  of  very  large  magnitude. 
This  is  particularly  critical  in  stations  with  very  short  records  like  the  North  Fork  of  the 
Virgin  River  below  Bullock  Canyon  (ID:  09405420),  for  which  the  graphical  mean  (200  cfs) 
represents  only  61%  of  the  arithmetic  mean  (326  cfs).    Conversely,  for  stations  with  long 
systematic  records  the  arithmetic  average  and  the  graphical  mean  tend  to  coincide. 

In  addition  to  the  MAPs,  graphical  estimates  of  flow  levels  for  recurrence  intervals  of 
5  yr  and  10  yr  were  also  obtained  from  the  single-site  plots.    Table  3  provides  a  summary  of 
all  the  information  graphically  derived.    Single-sites  frequency  analysis  were  conducted  for 
28  out  of  the  39  records.    Stations  with  very  short  records  were  not  considered.    The  left 
portion  of  Table  3  includes  basic  statistics  computed  from  the  systematic  records,  and  the 
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right  gives  flow  levels  for  three  different  quantiles.    The  significance  of  the  ratios  Q10/2.7  and 
Q5/2  7  will  be  discussed  latter  in  Section  3.5. 

Only  systematic  records  of  peak  flows  were  considered  for  the  frequency  analysis. 
This  study  recognizes  the  importance  of  historic  flows  in  the  computation  of  flood  frequency 
analysis.    For  instance,  Webb  and  Baker  (1987)  analyzed  the  effect  of  incorporating 
estimated  historic  floods  (floods  not  included  in  the  systematic  record)  in  the  computation  of 
the  flood  frequency  on  the  Escalante  River  (site  3  in  our  study).    However,  the  addition  of 
historic  peaks  (with  very  high  return  periods,  around  and  over  100  yr)  practically  doesn't 
change  the  distribution  of  peaks  for  low  return  periods  (in  the  1  to  5  yr  range),  as  it  can  be 
inferred  from  Webb  and  Baker's  study.    Furthermore,  a  basin-wide  effort  to  detect  and 
estimate  historical  peak  flows  is  beyond  the  scope  of  this  study. 

3.5       Langbein's  Homogeneity  Test 

A  customary  reference  to  analyze  hydrological  homogeneity  for  a  regional  flood 
frequency  analysis  is  the  test  developed  by  W.B.  Langbein  of  the  U.S.  Geological  Survey. 
Langbein  conceived  an  statistical  test  based  on  comparing  the  slope  of  frequency  curves  from 
individual  streams  within  a  study  region.    The  underlying  assumption  is  that  if  individual 
streams  having  a  wide  range  of  drainage  areas  have  frequency  curves  of  about  equal 
steepness,  the  records  of  peaks  will  be  homogeneous  with  respect  to  the  flood-producing 
mechanism.    The  variation  of  these  slopes  can  be  tested  to  see  whether  the  spread  could 
occur  by  chance  among  samples  from  the  same  population  due  simply  to  sample  variability. 
The  problem  is  then  to  calculate  the  spread  to  be  expected,  and  to  set  limits  to  the  spread 
that  is  acceptable.    Details  of  the  formulation  of  Langbein's  test  can  be  found  in  Dalrymple 
(1960). 

In  Langbein's  test  the  slope  of  the  frequency  curves  is  measured  by  the  ratio  of  the  10 
yr  to  the  2.33  yr  flood  (the  Gumbel  distribution  is  used  by  Langbein  to  represent  extreme 
events).    In  our  study  mean  annual  peaks  are  associated  with  a  slightly  different  recurrence 
interval,  2.7  yr.    The  acceptable  spread  of  the  ratio  is  associated  with  the  standard  deviation 
of  the  reduce  variate  of  the  Gumbel  distribution.    Typically  2  standard  deviations  are 
adopted.    This  implies  that  the  permissible  range  should  contain  95%  of  the  possible 
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observations.    The  10  yr  flood  has  been  used  to  compute  the  slope  since  this  is  the  longest 
recurrence  interval  for  which  many  records  will  give  dependable  estimates.    Flood  ratios 
have  been  computed  for  27  out  of  the  39  gaging  sites  based  on  the  extended  series  of  peak 
flows.    The  remaining  12  stations  have  records  too  short  to  be  considered  for  this  analysis. 
Sequential  steps  of  the  computational  procedure  are  included  in  Table  4.    The  last  two 
columns  in  Table  4  provide  the  plotting  coordinates  for  the  stations  in  Langbein's  test  chart. 
Tr  is  the  recurrence  interval  estimated  from  the  single-site  frequency  curve  for  a  flow  level 
given  by  the  product  of  Q2.7  times  the  average  ratio.    The  adjusted  period  of  record  is 
computed  as  the  length  of  the  systematic  record  plus  half  the  length  of  the  extended  portion 
of  the  record. 


Table  4.    Data  for  Langbein's  Homogeneity  Test 


Site 
No. 


USGS 
I.D. 


Mean  Ann. 
Q2.7 


Ratio 
Q10/Q2.7 


Q2.7 
x  Avg.  Rat 


Tr  for  prev. 
column 


Adjusted 
Per.  Rec. 


1 

09336400 

1335 

2  30 

2739 

8.0 

26 

2 

09337000 

240 

2.25 

492 

8.5 

45 

3 

09337500 

1430 

1.91 

2934 

12.9 

45 

4 

09338900 

910 

2.41 

1867 

7.3 

26 

6 

09379800 

1885 

2.19 

3867 

8.7 

14 

7 

09381100 

1165 

2.26 

2390 

8.3 

30 

8 

09381500 

3550 

1.75 

7284 

14.0 

32 

9 

09381 800 

2925 

2.12 

6001 

9.1 

29 

10 

09403500 

1290 

1.53 

2647 

60.0 

26 

11 

09403600 

850 

2.31 

1744 

78 

30 

12 

09403700 

1485 

1.51 

3047 

80.0 

26 

13 

09404450 

195 

2.41 

400 

7.2 

31 

14 

09404500 

1105 

2.18 

2267 

8.3 

25 

23 

09405500 

2285 

1.98 

4688 

11.5 

77 

24 

09405900 

2520 

1.96 

5170 

12.0 

22 

25 

09406000 

5075 

2.06 

10412 

10.0 

81 

27 

09406300 

220 

2.61 

451 

68 

29 

29 

09406700 

290 

2.52 

595 

7.1 

26 

31 

09408000 

655 

319 

1344 

5.4 

33 

32 

09408400 

95 

2.68 

195 

6.9 

35 

33 

09409500 

400 

2  70 

821 

6.5 

27 

34 

10174500 

720 

1  61 

1477 

50.0 

69 

35 

10241400 

100 

2.50 

205 

6.4 

28 

36 

10241470 

110 

2.41 

226 

7.3 

29 

37 

10241600 

160 

3.09 

328 

5.5 

29 

38 

10242000 

1150 

2.24 

2359 

8.2 

63 

39 

10241000 

430 

1.81 

882 

14.2 

26 

Average  Ratio 


2.21 
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The  test  chart  in  Figure  8  displays  the  set  of  test  curves  that  show  within  what  range 
of  recurrence  intervals  an  estimate  of  the  10  yr  flood  should  be  for  a  specified  length  of 
record.    Stations  that  plot  within  the  two  bounds  can  be  considered  to  be  homogeneous 
according  to  Langbein's  criterion,  and  therefore  can  be  considered  as  potential  candidates  for 
regional  flood  frequency  analysis.    Only  three  stations,  sites  10,  12  and  34  lay  outside  the 
confidence  bands.    The  three  stations  have  flood  ratios  noticeably  lower  than  the  average 
ratio  (2.21)  for  the  entire  set.    This  is  also  reflected  by  the  marked  downward  curvature  of 
their  frequency  curves  presented  in  Appendix  IV.    A  close  examination  of  the  annual  series 
of  flow  peaks  at  site  34,  Sevier  River  at  Hatch,  help  us  explaining  the  large  departure  in 
plotting  position  of  this  site  from  the  rest  of  the  stations.    It  was  observed  that  74%  of  the  69 
registered  peaks  occurred  from  late  March  to  early  June,  the  snow  melt  season  (most  of  them 
during  May).    This  finding  certainly  contradicts  the  initial  supposition  of  rainfall  being  the 
dominant  mechanism  producing  floods.    The  other  two  stations,  sites  10  and  12,  showed  a 
more  balanced  partition  between  snow  melt  and  rainfall  peaks,  however,  the  short  records 
available  at  both  sites  prevent  us  from  further  examination. 


100 


10      15      20     25      30      35     40     45      50      55     60     65      70      75      80      85      90 
Effective  Length  of  Record,  in  Years 

Fig. 8   Langbein's  Homogeneity  Test  Chart 
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3.6       Regional  Prediction  of  Mean  Annual  Peak 

Abundant  literature  can  be  found  since  the  early  1960's  that  analyzes  the  statistical 
relation  between  mean  annual  floods  estimated  at  several  sites  in  an  homogeneous  region  and 
explanatory  hydrological  variables.    The  procedure  consists  on  regionalizing  the  peaks  data 
on  the  basis  of  physiographic  characteristics  of  river  catchment,  characteristics  of 
precipitation,  or  other  meteorological  characteristics  of  various  types.    This  approach  is  also 
part  of  the  Generalized  Least-Square  Method,  currently  being  used  by  the  USGS  for  regional 
flood  frequency  analysis,  where  different  return  period  quantiles  are  related  to  basin 
characteristics  via  multiple  regression.    Based  on  results  from  previous  investigations,  this 
study  has  identified  a  set  of  hydrologic  variables  that  may  affect  the  occurrence  of  floods  in 
the  region  where  Zion  National  Park  is  located.    An  important  requirement  for  the  selection 
of  the  independent  variables  was  that  they  have  to  be  readily  obtainable  from  topographic  and 
meteorological  maps.    Compilation  of  the  hydro-meteorological  characteristics  for  all  the 
watersheds  was  a  tedious  and  time  consuming  task  despite  the  use  of  the  GIS  developed  for 
the  study  region.    Several  maps  of  the  watersheds  generated  during  this  stage  of  the  work  are 
included  in  Appendix  V.    The  maps  show:  boundaries  of  the  watersheds,  general  drainage 
network  and  the  adopted  main  water-course  (to  measure  channel  length  and  slope),  location 
of  the  basin  centroid  and  flow  gaging  site,  and  isopluvials  of  precipitation  intensity.    Table  5 
contains  a  summary  of  the  compiled  geographical  and  meteorological  information  for  all 
watersheds  with  the  exception  of  site  30  (site  30  was  only  used  for  peaks  data  extension 
purposes). 

A  brief  description  of  the  hydrological  characteristics  listed  in  Table  5  follows: 

Gage  Datum        (D)       indicates  the  elevation  of  the  terrain  where  the  USGS  gaging  station 

is  or  was  located.    Data  obtained  from  USGS  records. 

Watershed  Centroid        indicates  the  location  of  the  centroid  of  the  watershed.    Latitude  and 

longitude  coordinates  expressed  in  degree-minute-second.    Data 
obtained  from  the  GIS. 

Drainage  Area     (A)       gives  the  areal  extent  of  the  watershed  contributing  to  peak 

discharge, expressed  in  square  miles.    Data  obtained  from  the  GIS. 

30 


Table  5.  Geographical  and  Meteorological  Characteristics  of  Watersheds 

IIC/-e  RaalnM  Clittnn  

Watershed  f**»»*™i** 

Drainage 
Area  (GIS) 

Mean  Basin 

Channel  Ch 
El.®  85%L 

Channel 
Slope 

50%-€hr 
Preelp.lnt. 

2yr-24hr 

Preelp.lnt. 

■  ■-—■■  uoutj  <jayniy  *j 

Site  No      ID. 

Datum 

Latltud 

Longltud 

Elevation 

El.®  10%L 

10-85%Len  Total  Lngth 

(t) 

(mi2) 

(ft) 

(ft) 

(ft) 

(mi) 

(mi) 

(ft/mi) 

(inches) 

(inches) 

1 

09336400 

6250 

37  41  52 

1 1 1  48  34 

54  79 

7620 

6400 

7350 

9  48 

12  64 

100  2 

1  07 

1  59 

2 

09337000 

6400 

37  59:34 

111  39  29 

66  43 

8890 

6420 

9960 

13.60 

18  14 

260  2 

1  17 

1  95 

3 

09337500 

5670 

37  50  48 

111  44  21 

317  44 

8030 

5780 

7210 

17  24 

22  99 

82  9 

1  07 

1  76 

4 

09338900 

5700 

37  57  50 

1 1 1  22  52 

65  79 

7680 

5800 

8920 

12  59 

16  79 

247  8 

1  05 

1  70 

5 

09339200 

5080 

37  32  39 

111  31  55 

145.55 

6170 

5160 

6670 

24  38 

32  51 

61  9 

0  90 

1  39 

6 

09379800 

4240 

37  10  55 

1 1 1  49  40 

90  73 

5110 

4290 

5380 

16  36 

21  81 

66  6 

0  95 

1  30 

7 

09381100 

5980 

37  37  08 

1 1 1  52  56 

21  75 

7120 

6220 

7300 

11  27 

15  02 

95  9 

1  05 

1  52 

8 

09381500 

5440 

37  37  15 

112  00  14 

200  95 

6890 

5590 

7640 

19  32 

25  76 

106  1 

1  03 

1  46 

9 

09381800 

4300 

37  26:50 

1 1 1  59  40 

649  03 

6390 

4410 

6550 

48  35 

64  46 

44  3 

0  96 

1  38 

10 

09403500 

6220 

37  23  58 

11227  19 

73  51 

7250 

6240 

7480 

14  74 

19  65 

84  1 

1  19 

1  78 

11 

09403600 

5060 

37  15  53 

112  30  30 

192  64 

6670 

5130 

7060 

29  07 

38  76 

66  4 

1  12 

1  62 

12 

09403700 

5000 

37  12  39 

112  21  02 

245  92 

6300 

5200 

6560 

25  71 

34  28 

52  9 

1  05 

1  48 

13 

09404450 

5900 

37  24  50 

112  35  09 

74  30 

7300 

6030 

7140 

10  12 

13  49 

109  7 

1  17 

1  75 

14 

09404500 

5620 

37  15  23 

112  44  43 

5  57 

6110 

5610 

5950 

2  56 

3  42 

132  8 

1  12 

1  68 

15 

09404900 

3940 

37  15  59 

112  42  10 

343  48 

4200 

6760 

36  06 

48  08 

71  0 

1  11 

1  56 

16 

09405200 

7680 

37  33  34 

112  52  35 

6  84 

8484 

7780 

9080 

3  58 

4  78 

363  0 

1  42 

2  22 

17 

09405250 

7640 

37  31  52 

112  51  02 

9  33 

7818 

7720 

8650 

3  03 

4  03 

307  3 

1  45 

2  24 

18 

09405300 

8320 

37  31  38 

113  02  41 

10  22 

7592 

8420 

9580 

3  15 

421 

367  8 

1  32 

2  04 

19 

09405400 

7530 

37  29  46 

112  46  37 

5  65 

7650 

8560 

2  22 

2  95 

4106 

1  42 

2  22 

20 

09405420 

6420 

37  27  58 

11247  13 

29  51 

6590 

8150 

5  58 

7  44 

279  6 

1  37 

2  04 

21 

09405450 

6000 

37  26  56 

112  46  59 

45  35 

8010 

6120 

7940 

8  00 

10  67 

227  5 

1  35 

2  04 

22 

09405490 

4400 

37  26  17 

11254  09 

288  18 

4550 

7580 

19.35 

25  80 

156  6 

1  29 

1  97 

23 

09405500 

3970 

37  24  30 

112  54  11 

343  42 

7350 

4250 

7330 

24  47 

32  62 

125  9 

1  28 

1  94 

24 

09405900 

3680 

37:18:54 

113.05  41 

95  12 

6542 

3810 

6350 

12  15 

1621 

209  0 

1  12 

1  76 

25 

09406000 

3440 

37  18  10 

112  52  30 

956  40 

6400 

3580 

6610 

38  54 

51  39 

78  6 

1  24 

1  82 

26 

09406150 

3040 

37  22:52 

113:1001 

96  28 

3160 

6200 

23  80 

31  73 

127  7 

1  06 

1  66 

27 

09406300 

5680 

37  31  54 

11307  13 

9  84 

7950 

5880 

8040 

4  60 

6  13 

470  0 

1  12 

1  85 

28 

09406500 

4450 

37  29  34 

113  1541 

132  98 

4750 

7030 

1280 

17  07 

178  1 

0  99 

1  58 

29 

09406700 

5290 

37  22  54 

113  22  03 

11  13 

7210 

5430 

7200 

3  55 

4  74 

498  5 

1  22 

1  90 

30 

09406800 

4100 

6720 

31 

09408000 

4000 

37  18  42 

11324  51 

15  56 

6360 

4220 

6640 

6  35 

8  46 

381  2 

1  08 

1  69 

32 

09408400 

6640 

37  22  06 

113  26  35 

18  67 

8720 

6730 

9120 

4  73 

6  31 

505  3 

1  31 

2  14 

33 

09409500 

4800 

37  27  53 

113  48  35 

33  96 

6070 

5140 

6380 

9  29 

12  39 

133  4 

1  02 

1  57 

34 

10174500 

6870 

37  34  09 

112  36  41 

332  32 

8480 

6920 

8170 

22  07 

29  42 

56  6 

0  98 

1  54 

35 

10241400 

6740 

37  52  46 

112  39  12 

16  04 

7470 

6910 

8100 

5  33 

7  38 

223  1 

0  97 

1  59 

36 

10241470 

6900 

37  44  36 

112  47  42 

14  05 

8450 

6850 

10040 

5  93 

7  91 

537  9 

1  13 

1  87 

37 

10241600 

6313 

37  43  59 

112  53  56 

23  75 

8230 

6740 

9200 

7  08 

9  44 

347  5 

1  08 

1  58 

38 

10242000 

6000 

37  37  53 

11256  21 

80  39 

8640 

6200 

8290 

11  08 

14  77 

188  6 

1  14 

1  75 

39 

10242100 

5835 

37  35  26 

113  03  57 

12  94 

8032 

5910 

8030 

4  22 

5  62 

502  7 

1  19 

1  97 

Max 
Mm 


8320 
3040 


956  40 
557 


8890 
5110 


48  35 
2  22 


64  46 
2  95 


537  9 
44  3 


1  45 
0  90 


2  24 
1  30 


Mean  Basin  El.    (E) 


Channel  Length   (L) 


Channel  Slope     (S) 


measures  the  median  elevation  of  the  watershed,  in  feet.    It  can  be 
interpreted  as  an  index  of  the  orographic  effect  in  the  type  of 
precipitation  in  the  basin.    Data  obtained  from  the  USGS  and  GIS. 

measures  the  length  of  the  main  water  course  from  the  watershed 
outlet  to  the  basin  divide,  in  miles.    The  selected  channel  is  not 
necessarily  the  one  with  the  longest  path,  but  the  stream  suspected  to 
be  the  larger  contributor  of  peaks.    Data  obtained  from  the  GIS. 

indicates  the  average  rate  of  change  in  elevation  of  the  main  channel, 
in  feet  per  mile,  measured  from  the  10th  to  the  85th  percentile  of  the 
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total  channel  length  upstream  from  the  gage.    This  is  an  indirect 
measure  of  the  basin  gradient.    Data  obtained  from  the  GIS. 

Precip.  Inten.       (P6)      a  precipitation-frequency  index  based  exclusively  on  rainstorms  with 
(50%-6hr)  a  50%  exceedanceprobability  and  a  6-hour  duration.    The  value  of 

the  variable,  in  inches,  is  associated  with  a  point  probability  at  the 
basin  centroid.    Derived  from  NWS  map  using  the  GIS. 

Precip.  Inten.       (P24)     a  precipitation-frequency  index  representing  the  occurrence  of  all 
(2yr-24hr)  precipitation  events  with  a  recurrence  interval  of  2-years  and  24-hour 

duration.    P24  was  obtained  in  the  same  form  as  P6. 

Figure  9  shows  descriptive  plots  where  Q2.7,  the  explained  variable,  is  individually 
related  to  six  of  the  hydrological  factors  (explanatory  variables)  affecting  the  magnitude  of 
peaks.    They  are:  drainage  area,  gage  datum,  mean  basin  elevation,  channel  length,  channel 
slope  and  precipitation  intensity.    Overall,  the  six  bivariate  correlations  (in  the  logarithmic 
domain)  give  satisfactory  results.    In  particular,  log-Area,  log-Slope,  log-Length  and  Log- 
Datum  depict  strong  correlative  associations  with  log-Q2.7.    Mean  basin  elevation  and 
precipitation  intensity  show  somewhat  larger  scatter  than  the  previous  four  variables,  but  still 
clear  relationships  can  be  noted.    It  is  hoped  that  when  the  NWS  releases  the  new 
precipitation  maps  for  south-western  Utah  the  P24  relation  can  be  improved.    The  relation 
between  log-Q2.7  and  log-P6  (not  shown)  was  observed  to  be  the  weakest.    The  same  simple 
correlation  task  was  conducted  using  log-MAP  as  the  dependent  variable.    Results  were 
similar  to  those  found  for  Q2.7,  except  that  a  larger  number  of  sites,  38  instead  of  28,  were 
include  in  the  latter  analysis.    The  log-domain  is  used  for  the  plots  in  Figure  9  to  reflect  the 
same  type  of  relations  used  in  the  regional  multivariate  model  to  be  introduced  in  Section  4. 

It  should  be  pointed  out  that  the  basin  characteristics  mentioned  above  are  not 
independent  among  themselves.    Several  of  them  are  significantly  interrelated  such  as 
drainage  area  and  channel  length,  channel  slope  and  gage  datum,  etc.    Related  hydrologic 
characteristics  should  not  be  used  together  for  predicting  purposes.    Our  objective  now  is  to 
identify  the  most  efficient  combination  of  hydrological  factors  in  explaining  the  variations  in 
the  mean  of  peak  flows  within  the  study  region.    There  is  a  statistically  defined  limit  to  the 
number  of  variables  necessary  to  relate  the  mean  annual  peaks  to  basin  characteristics. 
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25  3  35 
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0  3  0  4  0  5  0  6  0  7 

Log(e)  Precip  Int  (2yr-24hr) 


Fig. 9  Descriptive  Plots  of  Mean  Annual  Peaks  (Q2.7)  vs.  Hydrological  Factors 
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The  independent  variables  are  selected  using  multiple  linear  regression  analysis  with  stepwise 
technique.    Variables  are  included  in  the  order  determined  by  the  amount  of  contribution  to 
variance  reduction.    The  F-value  is  the  criterion  used  to  decide  the  inclusion  of  a  hydrologic 
characteristic  in  the  regression  equation.    Variables  are  added  one  by  one  to  the  model  and 
the  F-statistic  for  a  variable  to  be  added  must  be  significant  at  a  predefined  level.    The 
stepwise  process  ends  when  none  of  the  variables  outside  the  model  has  a  significant  F- 
statistic.    Table  6  provides  a  summary  of  the  results  of  the  stepwise  procedure. 

Table  6.    Summary  of  Stepwise  Procedure  for  Dependent  Variable  MAP 


Step  1 

Variable  AREA  Entered 

R-square  =  0.63 

DF 

Sum  of  Squares 

Mean  Square     F 

Prob  >  F 

Regression 

1 

46.2338 

46.2338             60.47 

0.0001 

Error 

36 

27.5241 

0.76456 

Total 

37 

73.7579 

Step  2 

Variable  DATUM  Entered 

R-square  =  0.67 

DF 

Sum  of  Squares 

Mean  Square     F 

Prob  >  F 

Regression 

2 

49.6594 

24.8297             36.06 

0.0001 

Error 

35 

24.0985 

0.68853 

Total 

37 

73.7579 

Step  3 

Variable  PRE24  Entered 

R-square  =  0.70 

DF 

Sum  of  Squares 

Mean  Square     F 

Prob  >  F 

Regression 

3 

51.6040 

17.2013             26.40 

0.0001 

Error 

34 

22.1539 

0.65159 

Total 

37 

73.7579 

No  other  var 

iable  met  the  0.150  significance 

level  for  entry  into  the  model 

Notice  in  Table  6  that  70%  of  the  regional  variance  of  the  mean  of  log-peaks  is 
explained  by  the  first  three  factors  entered  into  the  model,  with  log-area  as  the  single 
dominating  factor.    All  variables  in  the  model  are  significant  at  the  0.150  level.    No  other 
variable  met  the  0.150  significance  level  for  entry  into  the  model.    This  finding  represents  a 
typical  result.    However,  A,  D  and  P24  should  not  be  interpreted  as  the  only  variables 
affecting  flow  peaks,  but  rather  they  represent  the  most  efficient  subset  of  basin  character- 
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istics  available  to  predict  peaks.    In  fact,  when  the  stepwise  procedure  was  repeated  to 
predict  Q2.7,  this  time  only  two  variables  were  entered  to  form  the  predictive  model,  S  and 
D,  explaining  71%  of  the  Q2.7  regional  variance.    Figure  10  is  another  way  to  represent  the 
results  of  the  multiple  regression  analysis  by  plotting  observed  versus  predicted  log-mean 
annual  peaks  (graphical  estimate  Q2.7).    Figure  10  shows  an  acceptable  scatter  of  the  points 
around  the  45°  line.    The  predicting  capability  of  the  model  can  still  be  improved  if 
limitations  are  imposed  on  the  range  of  those  basin  characteristics  considered  in  the 
regressions,  in  other  words,  if  a  more  homogeneous  region  is  defined.    In  summary,  four 
hydrologic  factors  were  identified  as  potential  predictors  of  mean  annual  peaks  in  the  Zion 
region,  they  are:  drainage  area,  gage  elevation,  channel  slope  and  precipitation  intensity. 
More  details  about  the  final  form  of  the  model  adopted  to  represent  mean-annual  peaks  are 
given  in  Sections  4.1  and  5.3. 


5  6  7  8 

Observed  Log  of  Q2.7 


10 


Fig.  10  Observed  vs.  Predicted  Log  of  Mean  Annual  Peaks  (Q2.7) 
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4.0  MULTIVARIATE  LOGNORMAL  MODELS 

4.1  Regional  Multivariate  Lognormal-2  Model 

Many  hydrologic  variables  are  measured  from  the  base  zero  and  have  no  upper 
bound.    This  causes  their  distributions  to  be  asymmetrical.    It  has  been  found  through 
extensive  experience  (Yevjevich,  1982)  that  the  lognormal  probability  distribution  fits  well 
these  positively  skewed  distributions  in  hydrology,  including  frequency  of  floods.    The 
lognormal  distribution  has  been  widely  used  in  its  univariate  form  for  flood  frequency 
modeling.    Specifically  within  the  study  area,  frequency  analysis  at  a  single  site  showed  the 
log-Pearson  type  III  and  lognormal  distributions  equally  suitable  for  representing  maximum 
annual  peaks  at  the  N.  Fk.  Virgin  River  near  Springdale  (ID:  09405500)  (Diaz,  1992).    The 
similar  fitting  of  these  two  distributions  should  be  expected  since  watersheds  in  the  proximity 
of  Zion  National  Park  are  located  in  a  region  of  the  country  where  the  generalized  skew 
coefficient  for  peak  flows  series  is  practically  zero  (U.S.  Water  Resources  Council,  1981). 

Although  alternative  distributions  have  been  found  in  some  cases  even  more 
appropriate  than  the  lognormal  distribution  for  flood  frequency  modelling,  they  are  very 
difficult  to  deal  with  in  their  multivariate  forms.    This  study  applies  the  outcome  of  a  recent 
research  project,  Pons  (1992),  who  developed  a  methodology  for  regional  flood  frequency 
analysis  based  on  multivariate  lognormal  models.    Pons  chose  the  multivariate  lognormal 
distribution  with  2  and  3  parameters,  MLN2  and  MLN3  respectively,  because  they  have  a 
relatively  simple  form  while  still  maintaining  the  most  important  characteristics  needed  for 
flood  frequency  analysis.    The  following  sections  present  excerpts  from  Pons'  dissertation, 
including  the  models  formulation  and  parameter  estimation  procedures.    The  interested  reader 
should  consult  the  aforementioned  reference  for  a  thorough  description  and  discussion  of  the 
definition  and  testing  of  the  models. 

4.1.1    MLN2  -  Model  Formulation 

Basically,  the  lognormal  distribution  represents  the  normal  distribution  of  the 
logarithms  of  a  selected  random  variable,  in  this  case,  annual  peaks.    That  is,  if  the 
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logarithms  of  annual  peaks  y=ln  (x)  are  normally  distributed,  then  the  variable  x  is 
considered  to  be  logarithmic-normally  distributed.    Distribution  theory  provides  the  following 
properties  of  the  marginals  of  the  multivariate  lognormal-2  distribution, 

Y.  =  ln(X.)  -  N  (]>p),afy))  (4.1) 

n,<»40>)  (4.2) 

o](x)  -,*>»**»  (^»-l)  (43) 

r\j(x)  =  eJ     -1  (44) 

Yj(x)  =  r\](x)+3r\j(x)  (4  5) 

where  \ij(x)   is  the  mean  at  site  j,  o?(x)   is  the  variance  at  site  j,  "n^  (x)    is  the 
coefficient  of  variation  at  site  j  and  y7  (*)   is  the  coefficient  of  skewness  at  site  j  of  the 
untransformed  flood  variable.    Likewise,  \ij  (y)   is  the  mean  of  the  log-transformed  variable 
at  site  j  and  o){y)   is  the  variance  of  the  log-transformed  variable.    From  Eqs.(4.4)  and 
(4.5)  we  can  see  that  if  two  marginals  at  different  sites  have  equal  values  of  the  parameter 
a)  (y)  ,  they  will  have  equal  coefficient  of  variation  and  equal  coefficient  of  skewness.    This 
constitutes  a  restriction  imposed  by  this  distribution,  characterized  by  a  fixed  relationship 
between  the  coefficient  of  variation  and  the  coefficient  of  skewness. 

Similar  to  other  existing  methodologies  dealing  with  regional  floods  analysis,  the 
formulation  of  this  model  assumes  that  the  mean  of  maximum  annual  peak  series  can  be 
related  to  geographic  and  climatological  characteristics  of  the  catchments.    The  decision  as  to 
which  basin  characteristics  to  consider  has  been  investigated  in  Section  3.6.    It  is  based  on 
the  contribution  of  each  added  parameter  to  the  reduction  of  error  variance  in  the  regression 
equations  for  the  mean  and  the  basin  characteristics  .    Pons  relates  the  mean  of  peak  flows  to 
basin  characteristics  by  building  a  regression  model  for  the  mean  \ij  (y)  of  log-transformed 
peak  flows  at  site  j  in  the  following  manner, 
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\ij(x)  =  eJ      2J      =  CAfSfPf  (4.6) 

where  three  of  the  four  hydrologic  factors  identified  in  Section  3.6  were  chosen:  Aj  is  the 
drainage  area  at  site  j,  Sj  is  the  main  channel  slope  at  site  j  and  p7  is  the  precipitation 
intensity  at  site  j.    It  should  be  noted  that  the  addition  of  extra  physical  characteristics  to  the 
model  translates  into  additional  parameters.    Adopting  the  assumption  that  o)(y)   is  constant 
for  the  region  and  equal  to  o2  (y)  ,  Eq.(4.6)  can  be  written  as 

\i.(y)  =  c  +  a  \n(Aj)  +  s  ln(S.)  +  p  ln(Py)  (4.7) 

As  indicated  earlier,  the  assumption  of  a  constant  value  of  the  parameter  a)  (y)   for 
the  region  implies  that  also  the  coefficient  of  variation  and  the  coefficient  of  skewness  should 
be  constant  throughout  the  region.    In  addition,  in  a  hydrological  region  where  the  flood 
producing  mechanism  is  basically  unique,  a  relationship  between  the  cross  correlation 
coefficients  of  the  logs  of  peak  flows  data  and  the  inter-gaging  station  distances  can  be  built. 
Pons  (1992)  uses  a  negative  exponential  function  of  the  form, 

where  b  is  a  regional  parameter  and  d±j  is  the  distance  between  stations  i  and  j .    Though 
not  included  here,  Pons  demonstrates  that  the  cross  correlation  between  the  transformed 
series  and  the  cross  correlation  between  the  untransformed  series  are  functionally  related  with 
the  regional  variance  of  the  logs  a2  (y)  .    All  the  above  assumptions  permitted  Pons  justify 
the  following  form  of  the  lognormal-2  multivariate  model, 


-|      4    -i(x-u.(y))'v-1(x-uiy))  (4  9) 

M  -  (2*)  2  \V\  2  e  2  V 


where  v  and  ji(y)  are  defined  below,  and  '  denotes  the  transpose  of  the  vector, 
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c  +  a  ln^)  +  s  \n(S{)  +  p  ln^) 
c  +  a  \n(A2)  +  5  ln(52)  +  p  ln(P2) 
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o2(y) 

where  m  is  the  number  of  gaging  stations.    In  conclusion,  the  untransformed  series  of  peak 
flows  will  be  multivariate  lognormal  with  marginal  distributions  having  the  mean  expressed 
as  a  function  of  basin  characteristics,  a  common  coefficient  of  variation,  and  a  common 
coefficient  of  skewness. 

4.1.2    MLN2  -  Preliminary  Estimation  of  Parameters 

The  parameters  of  the  model:  c,  a,  s,  p,  b  and  o2  (y)  ,  are  ultimately  estimated 
using  the  maximum  likelihood  procedure  presented  in  Section  4.1.3.    However,  since  the 
equations  derived  from  that  procedure  cannot  be  solved  explicitly,  the  optimization  algorithm 
set  up  to  search  for  the  optimum  set  of  parameters  needs  to  be  initialized  with  a  set  of  initial 
estimates  of  the  searched  parameters.    This  is  a  feasible  task  since  all  parameters  have  a 
statistical  meaning.    For  instance: 
>  initial  estimates  for  c,  a,  s  and  p  are  obtained  from  a  multiple  linear  regression  of 

the  mean  of  log-transformed  peak  flows,  ji(y)  ,  as  the  dependent  variable,  versus  the 

logarithm  of  area,  logarithm  of  slope,  and  logarithm  of  precipitation  intensity  as  the 

independent  variables. 
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parameter  b  is  estimated  as  the  slope  of  the  regression  line  of  the  cross  correlation 
coefficients  of  the  log-transformed  peaks  versus  inter-station  distances,  forcing  the 
intercept  through  zero  to  equal  one. 

the  variance  of  the  log-transformed  peak  flows  o2  (y)   is  initially  estimated  as  the 
weighted  average  of  the  individual  variances  of  the  log  transformed  data,  with  weights 
equal  to  the  lengths  of  the  individual  records. 


4.1.3    MLN2  -  Maximum  Likelihood  Estimation  of  Parameters 

Maximum  likelihood  estimation  produces  estimates  that  posses  desirable 
characteristics  such  as  consistency,  and  also  of  efficiency  and  asymptotic  normality  under 
certain  regularity  conditions.    For  the  MLN2  model  these  regularity  conditions  are  satisfied, 
allowing  the  derivation  of  the  asymptotic  distribution  of  the  estimators  based  on  the  density 
function  given  by  Eq.(4.9).    For  the  most  general  case  of  unequal  samples  size  and  when  the 
observations  are  multivariate  in  space  and  independent  in  time,  the  likelihood  function  is 
obtained  as  the  product  of  the  likelihood  of  the  individual  observations  (annual  peaks 
recorded  at  multiple  stations).    After  taking  logarithms  of  the  likelihood  function  and 
factorizing  o2  (y)    out  of  the  variance  covariance  matrix,  Pons  gives  the  following  log- 
likelihood  LL  function, 


»  w.ln(27i)      i 
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and: 


mt  number  of  stations  observed  in  year  t 

^(y)  vector  of  means  of  the  mt  stations  present  in  year  t 

Et  matrix  of  cross  correlations  among  the  mt  stations  present  in  year  t 


The  method  of  maximum  likelihood  consist  in  estimating  the  parameters  c,  a,  s,  p,  b 
and  o2  (y)  such  that  LL  is  maximized.    As  indicated  earlier,  maximization  of  the  log- 
likelihood  function  in  Eq.(4.12)  can  not  be  solved  explicitly.    The  Powell  optimization 
algorithm  is  utilized  to  search  for  the  location  of  the  maximum  of  LL,  thus  yielding  the 
maximum  likelihood  estimates  of  the  parameters  of  the  model. 

4.2       Multivariate  Lognormal-3  Model 

As  indicated  in  Section  4.1.1,  the  fixed  relationship  between  the  coefficient  of 
variation  and  the  coefficient  of  skewness  for  MLN2  may  become  an  unacceptable  condition 
for  some  hydrologic  regions.    The  problem  can  be  circumvented  by  the  MLN3  distribution, 
which  offers  a  greater  flexibility  at  the  expense  of  the  complexity  and  difficulty  in  estimating 
parameters. 

4.2.1    MLN3  -  Model  Formulation 

When  the  lower  bound  of  the  empirical  distribution  of  peaks  is  not  zero,  it  is 
necessary  to  modify  the  lognormal  distribution  function  by  introducing  a  new  parameter,  the 
lower  boundary.    Typically,  the  lower  boundary  is  not  known  a  priori,  therefore,  the 
relations  for  parameters  given  for  MLN2  cannot  be  applied  directly.    From  distribution 
theory,  the  properties  of  a  marginal  of  the  MLN3  distribution  are  known  as, 


Tj  -  MW  -  N^jfyhojiy))  (4M) 

My)+-r«/(y) 
\ij(x)  =  Xj+e        2  (4.15) 
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2  2  (4.16) 

T1,W  =— —  (4.17) 

1  + 1 

^a>)_3eo>)+2  4l8 

Yjfr)  =  ; — 

(e^-l)2 

where  a){x)  ,  ^-(x)   and  y-j(x)  are  as  defined  in  Section  4.1.1,  and  \jlj(x)    denotes  the 
mean  at  site  j  of  the  untransformed  series.    The  additional  parameter  per  marginal 
distribution,  denoted  by  t7,  is  shown  in  the  set  of  equations  above.    This  parameter  can  not 
be  regionalized  since  it  represents  the  lower  bound  of  each  marginal  distribution.    However, 
by  assuming  a  constant  variance  a2  (y)  and  a  constant  coefficient  of  variation  r\  (x) 
throughout  the  region,  Xj  can  be  estimated  by, 
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(4.19) 


While  the  mean  of  the  log-transformed  peaks  are  related  to  basin  characteristics  for 
MLN2,  in  MLN3  instead,  the  mean  of  peak  flows  in  the  original  domain  is  utilized.    The 
same  form  of  model  shown  in  Eq.(4.6)  is  used.    Furthermore,  the  lag-zero  cross  correlation 
of  the  untransformed  peak  flows  can  also  be  represented  by  the  same  relationship  used  for 
the  previous  model,  Eq.(4.8).    Under  all  these  assumptions,  Pons  proposes  the  following 
form  of  the  MLN3  model, 

Y  =  ln(X-t)  ~  MVN(^(y),V)  (4.20) 
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where  m  is  the  number  of  stations.    In  conclusion,  the  untransformed  series  of  peak  flows 
will  be  multivariate  lognormal-3  with  marginal  distributions  having  the  mean  expressed  as  a 
function  of  basin  characteristics,  a  common  coefficient  of  variation,  and  a  common 
coefficient  of  skewness. 

4.2.2    MLN3  -  Preliminary  Estimation  of  Parameters 

The  parameters  of  the  model:  c,  a,  s,  p,  b,  a2  (y)  and  r\  (x)  are  estimated  using  a 
pseudo-maximum  likelihood  procedure.    Similar  to  the  procedure  implemented  for  the  MLN2 
model,  initial  estimates  of  the  parameters  are  necessary  prior  to  the  implementation  of  the 
pseudo-maximum  likelihood  estimator.    Parameters  are  initially  computed  as  follows: 

>  c,  a,  s  and  p  are  obtained  as  for  MLN2,  except  that  in  this  case  the  log  of  the  mean 
annual  peaks  at  each  station  are  used. 

>  parameter  b  is  estimated  from  the  cross  correlation  coefficients  of  the  untransformed 
observations.    Since  Xj  is  initially  unknown,  data  transformation  is  not  possible. 

>  the  variance  of  the  transformed  peaks  is  obtained  from  a  weighted  mean  of  sample 
coefficients  of  skewness.    After  obtaining  the  weighted  coefficient  of  skewness, 
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a2  (y)    is  computed  from  Eq.(4.18). 
>  the  coefficient  of  variation    t\(x)    is  obtained  from  a  weighted  average  of  the  sample 

coefficients  of  variation  of  the  untransformed  peaks. 

4.2.3    MLN3  -  Pseudo  Maximum  Likelihood  Estimation  of  Parameters 

In  the  lognormal-3  distribution  each  location  parameter  t7  has  to  be  smaller  than  the 
lowest  observed  value  at  the  corresponding  site  j .    In  other  words,  the  range  of  the 
distribution  is  constrained  by  one  of  the  parameters.    This  violates  one  of  the  regularity 
conditions  needed  to  obtain  the  asymptotic  properties  of  maximum  likelihood  estimators. 
Most  of  the  asymptotic  properties  of  the  MLN2  are  lost  since  a  full  maximum  likelihood 
estimation  procedure  for  the  MLN3  is  not  feasible.    Nevertheless,  the  method  still  consists  of 
finding  the  parameter  set  that  maximizes  the  log-likelihood  function,  though  this  time 
constrained  on  one  of  the  parameters.    From  there  comes  the  name  of  pseudo  maximum 
likelihood  estimator.    Following  same  previous  notation,  the  density  function  of  the  MLN3 
model  is  given  by, 


M  =  (2«)  2  \V\  2  e  2  n^-T  (4>23) 


where  ln(x-i)   is  the  vector  of  elements  ln(xJ.-tJ)  ,  and  the  other  terms  are  as  previously 
defined.    For  the  case  of  unequal  sample  sizes,  Pons  (1992)  provides  the  log-likelihood 
function  for  the  MLN3  model, 
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where  mt,  ^.(y)   and  Et  are  as  defined  below  Eq.(4.13).    Finally,  the  method  of  maximum 

likelihood  consists  in  estimating  the  parameters  c,  a,  s,  p,  b,  a2  (y)  and  n.  (x)   such  that  ll 
is  maximized.    As  indicated  earlier,  the  location  parameters  of  the  marginal  distributions 
need  to  be  constrained  so  the  maximum  values  they  can  take  correspond  to  the  minimum 
observed  values  at  the  gaging  station.  Since  lower  bounds  x7  are  not  explicitly  defined  as 
parameters  of  the  model,  they  become  constraints  to  the  rest  of  the  parameters  during  the 
optimization  process.    Again  the  Powell  algorithm,  but  this  time  with  some  modifications  to 
accommodate  for  the  constrained  system,  is  implemented  to  search  for  the  maximum  of  the 
log-likelihood  function. 
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5.0  MODELS  APPLICATION 

5.1  Validation  of  Underlying  Assumptions 

The  formulation  of  the  MLN2  and  MLN3  regional  models  is  based  on  two  major 
assumptions:  1)  the  variance  of  the  log-transformed  peaks  a2  (y)  is  constant  throughout  the 
study  region,  which  in  turn  implies  a  constant  coefficient  of  variation  n,  (x)   and  a  constant 
coefficient  of  skewness  y  (x)   for  the  region,  and   2)  the  regional  cross-correlation  (space 
correlation)  of  maximum  annual  peaks  between  pairs  of  stations  can  be  represented  by  an 
exponentially  decaying  function.    All  these  assumptions  are  being  graphically  examined  in  the 
following  sections. 
5.1.1    Regional  Characteristics  of  the  Coefficient  of  Variation 

Several  studies  from  the  early  1960's  to  the  present  have  assumed  the  coefficient  of 
variation  of  annual  peak  series  to  be  a  constant  within  a  region.    On  the  other  hand,  Benson 
(1962)  argues  that  the  coefficient  of  variation  depends  on  basin  characteristics  such  as 
drainage  area,  channel  slope  and  many  others.    Pons  (1992),  in  line  with  the  outcomes  from 
the  latest  studies,  avoided  an  explicit  formulation  of  the  coefficient  of  variation  in  his 
multivariate  models  (see  Section  4.1.1  and  4.2.1.).    In  order  to  comply  with  Pons' 
assumptions  it  is  necessary  that  we  test  the  hypothesis  indicating  the  relationship  between  the 
coefficient  of  variation  and  hydrological  characteristics  of  the  watersheds  as  non-existent  or 
at  least  negligible  within  the  region.    No  analytical  procedure  is  introduced  for  this  purpose, 
only  judgment  based  on  graphical  observations  is  presented. 

Figure  1 1  shows  four  descriptive  plots  of  the  coefficient  of  variation  of  the 
untransformed  peak  series  (from  Table  3)  versus:  a)  drainage  area,  b)  stream  gage  elevation, 
c)  channel  slope  and  d)  precipitation  intensity  (from  Table  5).    This  is  the  same  set  of 
hydrologic  variables  that  yielded  the  most  efficient  correlative  association  with  the  mean 
annual  peaks  (see  Section  3.6).    Note  in  Figure  11  that  Cv  values  in  all  four  plots  are  widely 
scattered  over  the  whole  range  of  the  hydrological  variables.    The  visual  analysis  indicates 
that  the  coefficient  of  variation  is  independent  of  the  analyzed  hydrological  characteristics, 
and  therefore  Cv  can  be  assumed  approximately  constant  within  the  region.    Any  small 
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Fig.  1 1    Coefficient  of  Variation  versus  Hydrological  Characteristics 
a)  Drainage  Area,   b)  Gage  Datum,    c)  Channel  Slope,   d)  Precipitation  Intensity 


dependance  of  the  coefficient  of  variation  on  basin  characteristics  such  as  for  channel  slope 
can  be  considered  negligible.    The  size  of  watersheds  included  in  the  analysis  can  be 
properly  limited  to  avoid  any  dependence  of  Cv  with  geographical  variables. 

5.1.2    Regional  Characteristics  of  the  Coefficient  of  Skewness 

The  typical  large  variance  found  on  the  estimation  error  of  the  skewness  coefficient 
has  prompted  many  investigators  to  regionalize  this  high  order  statistical  property  of  the 
annual  series  of  peaks.    The  idea  of  assuming  a  regional  coefficient  of  skewness  has  also 
been  proposed  by  the  U.S.  Water  Resources  Council  in  their  log-skew  regional ization 
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procedure  presented  in  Bulletin  17B  (1981).    Bulletin  17B  suggests  regional ization  of  the 
skewness  coefficient  when  the  sample  length  is  below  a  certain  threshold  level.    The  majority 
of  the  stations  within  and  around  Zion  National  Park  have  short  records  that  would  require  a 
regionalization  of  the  skewness  coefficient.    For  the  purpose  of  regionalization  the  Bulletin 
provides  a  generalized  skew  coefficient  map  derived  from  multiple  long-term  records.    It 
should  be  pointed  out  that  Bulletin  17B  makes  no  distinction  between  annual  series  with  snow 
melt  or  rainfall  peaks  for  the  derivation  of  the  generalized  skew  coefficient.    The  importance 
of  this  consideration  within  the  Zion  region  has  been  discussed  in  Section  3.1. 

Coefficients  of  skewness  estimated  from  the  log-transformed  systematic  records  at 
several  gaging  sites  are  shown  in  Figure  12  (see  Log-skew  data  in  Table  3).    Two  clusters  of 
points  can  be  distinguished.    A  first  group  for  stations  1  to  14  and  a  second  group  for 
stations  23  to  39.    The  two  groups  of  points  display  different  degree  of  scatter  around  a 
constant  value  very  close  to  zero  (-0.05).    The  group  on  the  left,  with  the  larger  scatter, 
represents  stations  located  east  of  the  Park.    Stations  in  this  group  have  short  and  very  short 
systematic  records  which  is  reflected  by  their  larger  variance.    Contrarily,  the  stations  in  the 
second  group,  in  close  vicinity  to  the  Park,  show  a  smaller  variability  around  the  average 
value.    This  group  contains  the  four  longest  records  available  in  the  region,  sites  23,  25,  34, 
and  38.    A  third  group  of  stations,  sites  15  through  22,  have  not  been  plotted  for  having  very 
short  records..    These  are  stations  located  mostly  in  the  headwaters  of  the  North  Fork  of  the 
Virgin  River. 

It  is  important  to  emphasize  the  critical  effect  that  the  length  of  a  record  has  in  the 
estimation  of  a  log-skew  coefficient.     Notice  for  that  purpose  the  plotting  position  of  stations 
with  long  records,  sites  2,  23,  25,  34  and  38,  all  marked  with  encircled  points.    These  five 
stations  represent  a  wide  range  of  watershed  sizes,  from  15  to  956  mi2,  with  gages  located 
within  an  also  wide  range  of  elevations,  from   6,400  to  3,440  ft.    The  five  stations  are 
among  the  closest  points  to  the  constant  average,  suggesting  that  as  more  long  records 
become  available,  stations  than  now  depart  from  the  constant  average  will  tend  to  converge 
toward  the  center  of  the  plot.    In  summary,  it  is  concluded  that  long-term  skew  coefficients 
in  the  Zion  area  are  better  described  by  a  regional  constant  value  rather  than  by  at-site 
estimates. 
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Fig.  12    Log-skew  Coefficients  in  the  Zion  Region 


5.1.3    Regional  Characteristics  of  the  Cross-Correlation 

The  relation  of  the  cross-correlation  coefficient  of  annual  hydrologic  series  to  inter- 
station  distance  has  been  extensively  studied  by  other  investigators,  see  for  instance 
Berndtsson  (1987).    Meteorological  events  that  have  a  noticeable  impact  on  hydrological 
series  are  regional  in  nature.    For  instance,  extreme  precipitation  events  and  their 
corresponding  peak  flows  is  a  typical  example  of  this  type  of  phenomenon.    Sometimes 
weather  systems  moving  over  a  region  create  response  patterns  that,  because  of  the  vast 
coverage  and  motion  of  the  system,  the  hydrological  response  vary  smoothly  over  the  area 
affected  by  the  event.    Several  examples  of  this  regional  conditions  can  be  identified  in  the 
Zion  region.    For  instance,  the  cloudburst  storm  that  caused  the  flood  of  December  6^  1966 
in  Southwestern  Utah  was  a  slow  moving  system  traveling  eastward  across  the  state.    Record 
24-hour  amounts  of  precipitation  were  measured  at  Alton,  Orderville  and  Kanab  (Butler, 
1972).    Peak  discharge  of  Virgin  River  at  Virgin  (site  25),  the  highest  of  record,  reached 
22,800  cfs.    Peak  discharges  were  also  the  highest  of  record  at  five  other  neighboring  gaging 
sites  as  a  result  of  the  same  storm,  sites  23,  25,  29,  30,  32  and  33.    Obviously,  under  this 
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conditions,  the  customary  assumption  of  independence  of  flood  events  is  not  satisfied.    This 
emphasizes  the  importance  of  utilizing  a  model  that  addresses  the  effect  of  interstation 
correlation  for  regional  flood  frequency  analysis. 

Estimation  of  the  cross-correlation  coefficient  is  carried  out  between  pairs  of  stations 
with  concurrent  records.    Unfortunately,  records  rarely  coincide  in  lengths  and/or  initiation 
dates,  hence  yielding  highly  variable  estimates  of  the  cross-correlation  coefficients.    Figure 
13  shows  intersite  dependance  of  peak  flows  derived  from  ten  sites  in  the  region.    The  sites 
are:  3,  11,  16,  17,  23,  24,  27,  28,  32  and  37.    They  were  selected  to  cover  a  wide  range  of 
interstation  distances.    As  expected,  the  cross-correlation  approaches  one  for  very  close 
stations,  and  decreases  exponentially  as  the  interstation  distance  increases.    The  exponent  'b' 
defines  the  rate  of  decay  of  the  exponential  function,  being  pre-estimated  by  the  model  using 
the  Least-square  method. 
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Fig.  13   Cross-correlation  vs.  Distance  Relationship  for  10  Selected  Stations 

The  level  of  cross-correlation  between  stations  in  a  region  defines  the  effectiveness  of 
a  multisite/regional  approach.    For  instance,  if  all  stations  in  a  region  were  perfectly 
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correlated,  Pj #J-  =  1  for  all  pairs  of  stations  i,  j  (and  of  equal  record  length),  there  would 
be  no  additional  information  to  be  extracted  from  neighboring  stations  to  increase  the 
reliability  of  a  quantile  estimate  at  a  given  site.    Conversely,  if  there  were  no  cross- 
correlation  of  peaks  among  the  sites,  the  effective  data  set  can  be  expanded  from  that  at  the 
single  site  to  the  entire  region.    In  practice  the  gain  due  to  regional ization  will  lie  somewhere 
between  no  gain  and  the  full  gain  possible  obtainable  only  when  there  is  no  cross-correlation 
and  a  strong  assumed  structure  among  site  parameters  (NRC,  1988). 

5.2       Selecting  an  Homogenous  Region  Based  on  Preserving  the  Model  Structure 

The  first  criterion  used  in  this  study  for  selecting  flow  gaging  sites  for  the  regional 
analysis  was  following  Thomas  and  Lindskow  (1983)  classification  of  flood  regions  for  the 
State  of  Utah,  see  Section  3.1.    Later  in  Section  3.5,    the  early  homogeneity  test  used  in 
conjunction  with  index  flood  methods  was  implemented.    Basically,  Langbeins'  homogeneity 
test  checks  for  similarly  shaped  flood  frequency  curves.    Based  on  Langbeins'  criterion,  only 
three  out  of  the  39  sites  originally  selected  from  Thomas'  Low  Plateau  region  were  found 
unfit  for  the  regional  analysis.    More  recent  regionalization  studies  have  focused  on  the  idea 
that  hydrologically  homogeneous  regions  need  to  be  delineated  in  terms  of  hydrologically 
meaningful  basin  characteristics  and  flood  statistics  as  well  as  in  terms  of  geographic 
locations,  see  for  instance  Wiltshire  (1986). 

This  section  proposes  an  alternative  approach  to  select  an  homogeneous  region.    It  is 
based  on  the  observed  series  of  peaks,  and  consist  in  analyzing  the  effect  that  the  series  have 
in  the  assumed  regional  structure  of  peaks  flows  as  defined  in  Sections  4.1.1  and  4.2.1.    In 
his  work  Pons  (1992)  analyzed  three  types  of  errors  that  can  cause  the  maximum  likelihood 
(ML)  procedure  to  yield  improper  values  of  the  model  parameters.    They  are:  1)  errors  in 
the  definition  of  the  mean  of  log-transformed  peaks  cr(m),    2)  errors  on  the  assumption  of 
constant  variance  of  log-peaks  <r(v),  and   3)  errors  on  the  assumption  of  a  fixed  relationship 
of  cross-correlation  versus  distance  cr(r).    After  comparing  these  three  types  of  errors,  Pons 
determined  that  deviations  from  the  assumed  relationship  between  the  mean  of  log-peaks  and 
basin  characteristics  have  the  most  detrimental  effect  on  the  maximum  likelihood  estimation 
procedure. 
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The  limitation  mentioned  above  can  be  circumvented  if  deviations  from  the  assumed 
regional  structure  of  peaks  are  minimized.    In  other  words,  we  should  find  a  group  of 
stations  that  satisfy  the  level  of  homogeneity  needed  for  the  ML  procedure  to  perform 
satisfactorily.    The  procedure  implemented  for  this  purpose  consists  in  running  the  model 
repeatedly,  dropping  one  station  at  a  time.    The  station  being  removed  from  the  set  at  each 
step  is  the  one  that  grants  the  largest  reduction  of  the  variance  of  the  error  term  under 
consideration.    Since  there  are  three  error  terms  to  minimize,  the  procedure  is  repeated 
entirely  for  each  one  of  them.    Figure  14  shows  the  results  of  the  series  of  runs.    The  three 
curves  display  the  error  terms  variance  as  a  function  of  the  number  of  stations  remaining  on 
the  set.    In  general  the  three  curves  display  a  smooth  decrease  of  the  error  terms.    Only  cr(v) 
decreases  more  rapidly  during  the  removal  of  the  first  six  stations. 


39      37      35      33      31      29      27      25      23      21      19      17      15      13      11       9        7 
Number  of  Stations  in  the  Set 

Fig.  14   Decrease  in  (r(m),  <r(r)  and  (r(v)  as  a  Function  of  Number  of  Stations 


It  should  be  mentioned  that  the  structure  of  the  cross-correlation  function  was  also 
investigated  using  the  interdistance  between  watersheds  centroids  instead  of  gaging  sites.    No 
improvement  in  the  relation  was  detected.    Intergage  distance  was  finally  adopted  since  it  is  a 
type  of  data  more  readily  available.    The  relation  of  the  mean  of  log-peaks  was  investigated 
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using  two  alternative  sets  of  hydrological  factors:  [A,  S,  P]  and  [A,  D,  P].    Again,  the 
results  showed  no  significant  differences.    The  set  [A,  D,  P]  was  finally  adopted  since  gage 
elevation  (datum)  is  easier  to  obtain  than  channel  slope. 

The  curves  in  Figure  14  show  how  the  variance  of  the  three  errors  decrease  as 
stations  are  progressively  removed  from  the  model.    However,  they  do  not  indicate 
specifically  which  stations  remain  in  the  set  at  each  level.    They  only  tell  us  that  for  instance, 
by  reducing  the  number  of  stations  from  38  to  15,  <r(m)  drops  from  0.84  to  0.14.    The 
additional  information  is  provided  by  Table  7  which  classifies  the  sites  based  on  their 
contribution  to  the  errors  variance.    A  station  displaying  a  dark  area  indicates  that  its 
presence  on  the  set  contributes  to  reduce  the  variance  of  the  respective  error.    Conversely,  a 
clear  area  indicates  that  the  presence  of  the  station  drastically  increases  the  error  term. 
Intermediate  cases  are  also  indicated. 

The  results  in  Table  7  were  used  to  select  the  most  desirable  subset  of  stations 
necessary  to  compose  the  regional  model.    The  relative  goodness  of  the  regional  structure  is 
judged  based  on  the  variance  of  the  three  errors  cr(m),  <r(r)  and  (^(v)  when  considered 
simultaneously.    Although,  as  indicated  earlier,  reducing    (^(m)    has  the  highest  priority. 
From  Table  7  we  can  conclude  that,  in  general,  the  model  rejected  first  the  stations  located 
the  farthest  from  the  Park,  as  indicated  by  the  larger  number  of  clear  areas  at  the  top  and 
bottom  of  the  table.    A  cluster  of  the  most  desirable  stations  (darker  areas  all  across)  is  found 
in  the  center  of  the  table,  from  sites  16  to  25.     The  stations  in  this  last  group  are  the  closest 
to  the  Park.    There  are  also  a  few  other  potential  good  sites  in  terms  of  reducing  the  errors 
variance,  like  sites  28  and  32. 

The  concept  of  areal  independence  was  another  criterion  for  selecting  the  final  subset 
of  stations  for  the  flood  regional  analysis.    Gaging  sites  are  considered  areally  independent 
when  areas  assigned  to  specific  gaging  sites  do  not  overlap,  i.e.,  individual  stations  represent 
mutually  exclusive  areas.    Some  of  the  stations  used  in  our  study  allow  for  partial 
overlapping  of  the  contributing  drainages.    See  for  instance  gages  at  sites  19,  20,  21,  22,  23, 
all  located  along  the  same  water  course.    Area  overlap  was  restricted  to  a  maximum  of  one 
third  of  the  larger  catchment.    In  summary,  based  on  the  error  analysis  described  above,  the 
length  of  records  available,  and  the  areal  independence  of  the  watersheds;  a  subset  with  ten 
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Table  7.    Classification  of  Gaging  Stations  Based  on  Their  Contribution  to  Errors  Variance 


Lowest  Error  Variance 
Intermediate  Error  Variance 
Highest  Error  Variance 
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gaging  stations  was  selected  for  estimating  the  parameters  of  the  regional  models  .  The  ten 
stations  are  marked  with  a  cross  in  Table  7.    They  cover  a  wide  range  of  drainage  areas, 
from  7  to  956  mi2,  as  well  as  mean  basin  elevations,  with  gaging  sites  ranging  from  3,500  to 
7,689  ft  in  elevation.    In  the  case  of  site  11,  although  it  exhibits  unfavorable  conditions  for 
reducing  cr(v),  it  was  still  included  in  the  subset  given  its  proximity  and  similar  topographic 
characteristics  with  site  15  (ungaged  site).    The  selected  subset  of  stations  have  specific 
common  properties  that  distinguish  the  subregion  from  the  neighboring  area.    In  other  words, 
from  the  point  of  view  of  the  structure  of  the  regional  model,  the  subregion  can  be 
considered  homogeneous.    The  reader  should  notice  that  the  selected  ten  stations  encircle 
several  potential  sites,  within  and  outside  the  Park,  where  frequencies  of  peak  flows  may 
need  to  be  estimated. 

5.3       Estimation  of  the  Parameters  of  the  Models 

The  maximum  likelihood  estimation  procedure,  briefly  described  in  Sections  4.1.3  and 
4.2.3,  was  applied  to  the  subset  of  stations  selected  in  the  previous  section.    Multiple  runs 
were  carried  out  to  explore  the  performance  and  sensitivity  of  the  ML  parameter  estimation 
procedure.    The  goodness  of  fit  of  the  MLN2  and  MLN3  models  was  judged  by  graphically 
comparing  the  resulting  regional  flood  frequency  curves  with  the  empirical  frequency 
distributions  at  the  stations.    Three  stations  were  selected  for  validating  the  regional  models: 
sites  13,  23  and  25.    They  were  chosen  due  to  the  good  quality  of  the  data,  extended  period 
of  record  (particularly  sites  23  and  25),  and  wide  range  of  watersheds  areas  and  elevations. 
Only  one  station,  site  13,  was  not  part  of  the  subset  of  stations  used  for  parameters 
estimation.    The  authors  recognize  that  the  validation  of  the  models  with  stations  not 
participating  in  the  parameter  estimation  procedure  would  have  been  preferable.    However, 
the  scarcity  of  stations  with  long  records  around  Zion  made  that  approach  virtually 
impossible.    Results  of  the  computations  are  presented  in  Table  8.    Table  8.  A  lists  eight  of 
the  selected  ten  stations  (and  their  corresponding  watershed  characteristics)  that  were 
ultimately  used  for  estimating  the  parameters  of  the  models  (sites  17  and  28  were  removed). 
It  also  shows  the  standard  errors  cr(m),  a2^)  and  <r(v)  as  estimated  initially  by  the  model. 
The  three  error  values  can  be  compared  with  the  curves  presented  in  Figure  14. 
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Table  8.    Summary  of  Results  from  ML  Procedure 
8.  A   Hydrological  Characteristics  of  Selected  Watersheds 


Stat    Map    USGS         Name 
No.    No.     ID 


Area     Datum     Ppt24 


Latitude      Longitude 


16 
19 

21 
23 
24 


9403600 
9405200 
9405400 
9405450 
9405500 
9405900 


Kanab  /Ka 
Deep  Cr 
NFVR  /Gle 
NFVR  /Zio 
NFVR  /Spr 
North  Cr 


7  25       9406000     Virgin  Rv 

8  32       9408400    Santa  Cla 


192.64 

6.84 

5.65 

45.35 

343.42 

95.12 

956.40 

18.67 


Sd.Dev  Residuals  for  Mean  of  Peaks  : 
Sd.Dev  Residuals  for  Cross  Correlation 
Sd.Dev  Residuals  for  Variance  of  Peaks 


5060. 
7680. 
7530. 
6000. 
3970. 
3680. 
3500. 
6640. 

0.358 
0.316 
0.293 


1.62 

2.22 
2.22 
2.04 
1.94 
1.76 
1.82 
2.14 


37.1003 
37.5217 
37.4717 
37.3903 
37.2100 
37.2369 
37.1981 
37.3839 


112.5483 
112.8833 
112.7794 
112.8250 
112.9778 
113.1506 
113.2075 
113.4819 


8.B    Estimated  Parameters,  ML  Procedure 


Param.  Description 


Lw.Bnd    Up.Bnd 


MLN2 

Initial  Final 


MLN3 

Initial  Final 


P(l)      Constant 

-6.00 

45.00 

25.849 

32.6451 

23.909 

23.9790 

P(2)      Coef  Ln(Area) 

-6.00 

12.00 

0.706 

0.4999 

0.647 

0.5755 

P(3)      Coef  Ln(Datum) 

-6.00 

12.00 

-2.855 

-3.5118 

-2.473 

-2.4326 

P(4)      Coef  Ln(Ppt24) 

-6.00 

12.00 

1.774 

1.3766 

0.755 

0.7476 

P(5)      Var.  LN(peaks) 

0.001 

6.00 

0.720 

0.6333 

0.324 

0.5808 

P(6)      Coef  X-Correla 

0.001 

0.10 

0.022 

0.0473 

0.030 

0.0405 

P(7)      Coef  of  Variat 

0.001 

10.00 

0.879 

0.9167 

No. Iterations  : 

28 

5 

No. Function  Evaluations  : 

4186 

896 

Value  Log-Likelihood  : 

253.5 

1598.5 

Table  8.B  provides  a  summary  of  the  computations  of  the  parameters  estimated  by  the  ML 
procedure  for  the  two  regional  models.    Figure  15  (top  graph)  shows  the  type  of  fitting 
achieved  with  MLN2  and  MLN3  at  site  23.    The  dots  represent  the  empirical  frequency 
distribution  of  the  measured  annual  peaks,  whereas  the  fitted  models  are  as  indicated  by  the 
legend  on  the  graph.    Both  regional  models,  MLN2  and  MLN3,  provide  very  good  and 
similar  fitting  of  the  available  data.    Also  the  single-site  fitting  of  the  data  using  a  univariate 
lognormal-2  model  is  provided  for  comparison.    The  regional  models  give  an  improved 
fitting  compared  to  the  univariate  model  for  recurrence  intervals  greater  than  10  years.    For 
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Fig.  15    Regional  and  Univariate  Fitting  of  Annual  Peaks 

a)   at  Site  23,  N.Fk.  Virgin  River  nr  Springdale,  b)  at   Site  25,  N.Fk.  Virgin  River  at  Virgin, 

c)  at   Site  13,  E.Fk.  Virgin  River  nr  Glendale 
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recurrences  smaller  than  7  years  the  fittings  by  the  univariate  and  multivariate  models  are 
almost  identical.    Similarly  good  results  were  obtained  when  evaluating  the  performance  of 
the  regional  model  for  sites  13  and  25,  see  graphs  15  b)  and  c),  where  the  fitting  by  the 
MLN2  model  is  shown.    The  extended  series  of  peaks  at  site  13  were  used  to  estimate  the 
plotting  positions.    It  should  be  mentioned  that  the  use  of  other  short  records  for  validation  of 
the  models  is  of  questionable  value.    Empirical  frequency  distributions  made  up  of  short 
records  typically  have  their  plotting  positions  highly  distorted. 

The  validation  process  presented  above  assures  the  good  performance  of  the  models 
for  sites  located  in  the  intermediate  and  lower  reaches  of  the  upper  Virgin  River  system. 
Unfortunately,  there  isn't  any  station  with  long  records  in  the  headwaters  upstream  from  the 
Park  (site  13  has  only  27  years  of  record).    The  lack  of  long  records  in  the  headwaters 
prevent  us  from  presenting  a  full  validation  of  the  models  including  the  high  elevation 
catchments.    Nevertheless,  it  is  expected  that  the  regional  character  of  the  models  plus  the 
inclusion  of  sites  with  very  high  mean  basin  elevation  during  the  calibration  stage  (sites  16, 
19  and  21),  despite  of  their  short  records,  will  yield  acceptable  results  even  for  watersheds 
located  at  very  high  elevation. 

5.4       Application  of  the  Models  to  Ungaged  Watersheds 

One  of  the  main  purposes  for  developing  a  regional  flood  frequency  model  is  to  be 
able  to  estimate  flood  frequency  curves  at  ungaged  sites.    As  demonstrated  in  the  previous 
two  sections,  the  parameters  of  the  model  should  be  derived  from  a  set  of  stations  adjacent  to 
the  ungaged  site  in  order  to  preserve  as  much  as  possible  the  underlying  structure  of  the 
probability  model.    Furthermore,  it  is  recommended  that  the  basin  characteristics  of  the 
ungaged  site  be  similar  to  the  watershed  characteristics  of  the  sites  used  for  estimating  the 
parameters  of  the  model.    The  procedure  for  applying  the  regional  models  at  ungaged  sites 
within  Zion  National  Park  and  vicinity  areas  is  as  follows: 

1.  Define  the  set  of  stations  to  be  used  for  estimating  the  parameters  of  the  model.    This 
task  has  been  completed  and  the  results  presented  in  Section  5.2. 

2.  Estimate  the  parameters  of  the  model  using  the  maximum  likelihood  procedure.    This 
task  has  also  been  completed  and  the  results  presented  in  Section  5.3. 
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3.  Determine  the  basin  characteristics  of  the  ungaged  watershed.  They  are:  drainage 
area  (in  square  miles),  terrain  elevation  of  the  selected  river  site  (in  feet),  and  the 
2yr-24hr  precipitation  intensity  (in  inches)  estimated  at  the  centroid  of  the  catchment. 

4.  Compute  the  mean  u7  (y)  ,  the  standard  deviation  o^y)  ,  and  lower  bound  x7  (the 
latter  only  for  MLN3)  of  the  marginal  distribution  j  based  on  the  estimated 
parameters  of  the  regional  model: 

-  for  MLN2:  \ijiy)  =  c  +  a  In  (A,.)  +s  IniSj)  +p  ln(PJ.) 

Oj(y)  =  vra2T>T 

-  for  MLN3:  [i,- U)  =  eIC+alaU^  +sln(5'}  ** **<*/" 


Oj{y)  =\/  a2  (y) 

T)  (X)  \ 


XJ 


1 
eo2(y)  _!)    2 


Vj(x) 


[ijiy)  =  LogeHijiyj-Xj]  -  1  o2  (y) 

5.         Compute  the  frequency  distribution  of  peak  flows  at  the  ungaged  site  as, 

-  for  MLN2:  YT  =  \ij{y)  +  t  Ojiy) 

-for  MLN3:  XT  =  Tj. +  e[^(y)  +t0^(y)] 

where  t  is  the  standard  normal  deviate  obtainable  from  tables  of  the  area  under  the 
normal  curve  or  using  a  polynomial  approximation.. 

Site  15,  in  the  East  Fork  of  the  Virgin  River  near  Springdale,  was  chosen  to 
demonstrate  the  capability  of  the  models  to  estimate  percentiles  of  the  frequency  distribution 
at  an  ungaged  location.    Site  15  can  be  considered  ungaged  since  only  two  years  of  record 
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are  available.    The  computations  indicated  in  points  3,  4  and  5  were  carried  out  using  a 
spreadsheet.    Results  of  the  spreadsheet  computations  are  presented  in  Table  9.    The 
spreadsheet  is  provided  with  this  report  to  facilitate  the  final  computations  to  the  user.    The 
estimated  frequency  distribution  at  site  15  is  shown  graphically  in  Figure  16.    The  two 
regional  models  provide  very  similar  estimations  of  the  frequency  distribution. 

It  is  interesting  to  compare  the  results  of  the  estimated  flood  frequency  curve  at  site 
15  with  one  of  the  two  annual  peaks  measured  at  the  site  so  far.    For  instance,  a  peak  of 
considerable  magnitude  was  observed  during  water-year  1993.    If  we  enter  the  frequency 
curve  in  Figure  16  with  the  value  of  the  maximum  instantaneous  discharge  for  that  water- 
year,  estimated  in  1,940  cfs,  we  can  determine  that  the  measured  peak  has  an  estimated 
frequency  of  occurrence  close  to  3-years.    Frequency  estimates  at  regular  recurrence 
intervals  can  be  read  directly  from  Table  9,  or  estimated  from  the  curves  in  Figure  16. 
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6.0     SUMMARY  AND  CONCLUSIONS 

The  main  outcome  of  this  study  was  the  implementation  of  a  regional  flood  frequency 
analysis  for  Zion  National  Park  and  vicinity  areas.    It  is  based  on  a  multivariate  lognormal 
model  that  explicitly  accounts  for  the  cross  correlation  among  the  series  of  annual  peaks  in 
the  region.    This  distinctive  characteristic  of  the  model  is  critical  in  the  Zion  region  where 
the  highest  peaks  have  been  recorded  during  the  same  year  at  several  stations.    This  finding 
precludes  the  use  of  traditional  Index-flood  methodologies  since  the  basic  requirement  of 
independence  of  events  is  violated.    Another  major  obstacle  for  the  analysis  is  the  limited 
database  of  annual  peaks  available.    This  restricts  the  applicability  of  another  technique,  the 
Generalized  Least  Square  method,  that  as  proved  by  other  investigators  (Pons,  1992), 
becomes  unreliable  when  the  number  of  available  stations  with  long  records  decreases. 

Two  regional  flood  frequency  models  with  the  statistical  characteristics  explained  in 
Section  4  were  implemented,  named  MLN2  and  MLN3.    They  will  allow  the  user  to  estimate 
flow  quantiles  (maximum  instantaneous  peaks)  at  ungaged  river  sites  based  on  readily 
available  geographical  and  meteorological  information.    The  models  can  also  be  used  to 
increase  the  reliability  of  estimated  flow  quantiles  at  sites  with  very  short  records  like  in  the 
case  of  headwater  catchments. 

This  study  used  Thomas  and  Lindskov  (1983)  classification  of  flood  regions  in  the 
State  of  Utah  to  preselect  gaging  sites  with  approximately  uniform  hydrologic  characteristics. 
Additional  analysis  and  testing,  including  Langbein's  homogeneity  test,  provided  further 
insight  into  the  hydrological  characteristics  of  the  annual  series  of  flow  peaks  in  the  region. 
It  was  found  that  peak  series  display  mixed-population  characteristics  resulting  from  snow 
melt  and  rainfall  events,  with  the  latter  as  the  predominant  mechanism.    The  ratios  of 
maximum  instantaneous  peaks  over  mean  daily  flows  for  the  same  day  (an  indicator  of  the 
degree  of  flashiness  of  the  stream)  are  low  and  very  constant  during  the  snow  melt  season, 
and  large  and  highly  variable  during  the  months  with  convective  thunderstorms.    High 
intensity  and  short  duration  summer  storms  give  rise  to  sharp  increases  and  fast  decay  of  the 
hydrographs,  with  peaks  reaching  up  to  about  30  to  60  times  the  mean-daily  flow  depending 
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on  the  mean  elevation  of  the  catchment.    On  the  other  hand,  during  the  snow  melt  season, 
usually  occurring  from  late  March  through  June,  instantaneous  peaks  rarely  exceed  mean- 
daily  flows  more  than  4  times.    The  probability  distribution  of  the  number  of  peaks  occurring 
at  different  times  of  the  year  was  also  examined.    A  distinctive  seasonal  pattern  was 
observed.    The  study  provides  a  methodology  for  the  estimation  of  the  probability  of 
occurrence  of  any  number  of  peaks  during  any  given  period  of  the  year. 

The  parameters  of  MLN2  and  MLN3  were  estimated  (using  the  maximum  likelihood 
procedure)  based  on  the  observed  series  of  peaks  from  a  subset  of  stations  that  are  known  to 
preserve  the  assumed  structure  of  the  regional  models.    The  selected  subset  of  stations  have 
specific  common  properties  that  distinguish  the  subregion  from  the  neighboring  area.    In 
other  words,  from  the  point  of  view  of  the  structure  of  the  regional  model,  the  subregion  can 
be  considered  homogeneous.    The  mean  of  peak  flows  in  the  region  can  be  in  great  part 
explained  by  basin  characteristics,  among  them:  basin  area,  channel  slope,  gage  elevation  and 
precipitation  intensity. 

Both  regional  models,  MLN2  and  MLN3,  showed  very  good  and  similar  fitting 
capabilities  when  compared  with  empirical  frequency  distributions  at  three  selected  stations. 
The  three  stations  are  representative  of  catchments  at  low,  intermediate  and  high  elevation 
areas  of  the  Upper  Virgin  River  Basin.    The  model  validation  process  assured  the  reliability 
of  the  regional  models  when  used  for  sites  located  at  intermediate  and  lower  reaches  of  the 
Upper  Virgin  River  system.    The  performance  of  the  models  at  high  elevation  sites  could  not 
be  thoroughly  verified  given  the  lack  of  long  records  for  headwaters  stations.    Nevertheless, 
we  expect  MLN2  and  MLN3  to  yield  reasonable  estimates  of  flood  quantiles  even  for 
watersheds  at  very  high  elevation.    A  spreadsheet  is  provided  with  this  report  to  facilitate  the 
final  computations  for  estimating  percentiles  of  the  frequency  distribution  curves  at  ungaged 
locations.    It  is  recommended  to  compute  and  compare  frequency  curves  obtained  from  both 
models,  MLN2  and  MLN3.    In  case  of  a  noticeable  difference  between  the  two  models  it  is 
suggested  to  select  the  curve  provided  by  MLN2  which  has  consistently  provided  more 
accurate  fittings. 
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Appendix  I  :    List  of  USGS  1:24000  (7.5  minute)  Topographic  Maps 


No. 

Map  Name 

USGS 

Code 

No. 

Map  Name 

USGS 

Code 

1 

Alton 

D4 

37112 

47 

Henrieville 

E8 

37111 

2 

Antimony 

A8 

38111 

48 

Hildale 

A8 

37112 

3 

Asay  Bench 

E5 

37112 

49 

Horse  Mountain 

D6 

37111 

4 

Bald  Knoll 

C4 

37112 

50 

Horse  Flat 

C7 

37111 

5 

Barker  Reservoir 

H7 

37111 

51 

Hurricane 

B3 

37113 

6 

Big  Lake 

A6 

38111 

52 

Jacobs  Reservoir 

A5 

38111 

7 

Boulder  Town 

H4 

37111 

53 

Johnson  Lakes 

A3 

37112 

8 

Brian  Head 

F7 

37112 

54 

Kanab 

A5 

37112 

9 

Bridger  Point 

A7 

37111 

55 

Kanarraville 

E2 

37113 

10 

Bryce  Point 

E2 

37112 

56 

King  Bench 

G3 

37111 

11 

Bryce  Canyon 

F2 

37112 

57 

Kolob  Reservoir 

Dl 

37112 

12 

Bull  Valley  Gorge 

Dl 

37112 

58 

Kolob  Arch 

D2 

37112 

13 

Butler  Valley 

D7 

37111 

59 

Little  Creek  Mountain 

A2 

37112 

14 

Calf  Creek 

G4 

37111 

60 

Long  Valley  Junction 

D5 

37112 

15 

Calico  Peak 

C8 

37111 

61 

Lower  Coyote  Spring 

B7 

37111 

16 

Canaan  Creek 

F6 

37111 

62 

Lower  Bowns  Reservoir 

A3 

38111 

17 

Canaan  Peak 

E7 

37111 

63 

Maple  Ridge 

D7 

37113 

18 

Cannonville 

El 

37112 

64 

Mount  Carmel 

B6 

37112 

19 

Carcass  Canyon 

E5 

37111 

65 

Navajo  Lake 

E7 

37112 

20 

Cedar  City 

Fl 

37113 

66 

Nephi  Point 

B2 

37112 

21 

Cedar  Mountain 

El 

37113 

67 

New  Harmony 

D3 

37113 

22 

Central  West 

D6 

37113 

68 

Nipple  Butte 

B6 

37111 

23 

Clear  Creek  Mountain 

C7 

37112 

69 

Orderville 

C6 

37112 

24 

Cogswell  Point 

D8 

37112 

70 

Page  Ranch 

E4 

37113 

25 

Collet  Top 

D4 

37111 

71 

Panquitch  Lake 

F6 

37112 

26 

Cottonwood  Mountain 

H6 

37112 

72 

Paragonah 

H7 

37112 

27 

Cutler  Point 

B4 

37112 

73 

Parowan 

G7 

37112 

28 

Dave  Canyon 

F5 

37111 

74 

Parowan  Gap 

H8 

37112 

29 

Death  Ridge 

E6 

37111 

75 

Petes  Cove 

D5 

37111 

30 

Deer  Spring  Point 

C2 

37112 

76 

Petrified  Hollow 

A2 

37112 

31 

Deer  Creek  Lake 

A4 

38111 

77 

Pine  Hollow  Canyon 

Al 

37112 

32 

Deer  Range  Point 

CI 

37112 

78 

Pine  Point 

B3 

37112 

33 

Eightmile  Pass 

Bl 

37112 

79 

Pine  Lake 

F8 

37111 

34 

Elephant  Butte 

A7 

37112 

80 

Pintura 

C3 

37113 

35 

Escalante 

G5 

37111 

81 

Podunk  Creek 

D3 

37112 

36 

Fivemile  Valley 

B8 

37111 

82 

Pollywog  Lake 

A7 

38111 

37 

Flanigan  Arch 

F8 

37112 

83 

Posy  Lake 

H6 

37111 

38 

George  Mountain 

E4 

37112 

84 

Rainbow  Point 

D2 

37112 

39 

Glendale 

C5 

37112 

85 

Roger  Peak 

H5 

37111 

40 

Grass  Lakes 

H8 

37111 

86 

Seep  Flat 

E4 

37111 

41 

Grass  Valley 

D4 

37113 

87 

Signal  Peak 

C4 

37113 

42 

Griffin  Point 

G7 

37111 

88 

Skutumpah  Creek 

C3 

37112 

43 

Hatch 

F4 

37112 

89 

Slickrock  Bench 

D8 

37111 

44 

Haycock  Mountain 

F5 

37112 

90 

Smith  Mesa 

C2 

37112 

45 

Hebron 

F7 

37113 

91 

Smithsonian  Butte 

Al 

37113 

46 

Henrie  Knolls 

E6 

37112 

92 

Springdale  West 

Bl 

37112 
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No. 

Map  Name 

USGS 

Code 

93 

Sringdale  East 

B8 

37112 

94 

Steep  Creek  Bench 

H3 

37111 

95 

Stoddard  Mountain 

E3 

37113 

96 

Straight  Canyon 

D7 

37112 

97 

Strawberry  Point 

D6 

37112 

98 

Summit 

G8 

37112 

99 

Sweetwater  Creek 

G8 

37111 

100 

Temple  of  Sinawava 

C8 

37112 

101 

The  Guardian  Angels 

CI 

37112 

102 

The  Barracks 

B7 

37112 

103 

Thompson  Pass 

A4 

37112 

No.        Map  Name 


USGS   Code 


104 

Tropic  Canyon 

Fl 

37112 

105 

Tropic  Reservoir 

E3 

37112 

106 

Upper  Valley 

F7 

37111 

107 

Virgin 

B2 

37112 

108 

Webster  Flat 

E8 

37112 

109 

West  Clark  Bench 

A8 

37111 

110 

White  Tower 

B5 

37112 

111 

Wide  Hollow  Reservoir 

G6 

37111 

112 

Wilson  Peak 

F3 

37112 

113 

Yellowjacket  Canyon 

A6 

37112 

No. 

Additional  Maps 

USGS 

Code 

114 

Fivemile  Ridge 

G5 

37112 

115 

Little  Creek  Peak 

H5 

37112 

116 

Red  Creek  Reservoir 

G6 

37112 

117 

Sunset  Flat 

E3 

37111 
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Watershed  with 

Outlet  at  Site  No. 

USGS  7.5'  Topographic  Maps 

1 

16,  42,  79,  106 

2 

6,  35,  52,  83,  85,  111 

3 

<1>,  <2>,  5,  40,  99 

4 

7,  31,  62,  56,  94 

5 

17,  19,  25,  28,  29,  75,  86,  117 

6 

9,  15,  36,  50,  61,  68 

7 

17,  79,  47,  106 

8 

<7>,  10,  11,  12,  18,  89,  104 

9 

<7>,  <8>,    9,  13,  15,  32,  33,  36,  50,  61,  77,  81,  84,105,  109 

10 

1,  4,  38,  39,  60,  81 

11 

<10>,  27,  54,  64,  110,  113 

12 

1,  4,  27,  53,  66,  76,  78,  81,  88,  103 

13 

1,  39,60,  69,  97 

14 

23,  64,  69,  102 

15 

<13>,  <14>,  34,  93,  96,  110,  113 

16 

65,  108 

17 

<16> 

18 

21,  57 

19 

65,  96,  97 

20 

<19> 

21 

<19> 

22 

<16>,  <17>,  <19>,  23,  24,  69,  100,  101 

23 

<22>,  92,  93,  102 

24 

58,  90,  92,  101,  107 

25 

<13>,  <14>,  <23>,  <24>,  34,  48,  91,  93,  110,  113 

26 

51,  57,  58,  80,  90,  101,  107 

27 

21,  55 

28 

<27>,  41,  57,  58,  67,  70,  95 

29 

41,  67,  80,  87 

31 

80,  87 

32 

41,  87 

33 

22,  45,  63 

34 

1,  3,  8,  38,  43,  44,  46,  60,  65,  71,  96,  97,  105 

35 

26,  114,  115,  116 

36 

8,  73 

37 

8,  37,  73,  98 

38 

8,  20,  21,  37,  43,  108 

39 

21 

Note:    <x>  indicates  that  all  maps  from  watershed  V  are  also  included  in  the  current  watershed 
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Appendix  II:        Recorded  and  Extended  Series  of  Annual  Peak  Flows, 

Basic  Statitics 
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Gaging  Stations 


Upp. Valley       Pin*  Cr       EscaJarte       Deer  Cr      Tw»ntymil       Coyote       Henrlevllle    Paria  Can.     Parla  Kan.    Kanab  Qle    Kanab  Kan.     Johnson     EFVR  n.  01       MlneraJ       EFVR  Spr.      Deep  Cr     EF  Deep  Cr      Crystal       NFVR  Gie 
09338400  09337000  09337500  09338900  09339200  09379800  09381100  09381500  09381800  09403500  09403600  09403700  09404450  09404500  09404900  09405200  09405250  09405300  09405400 


Basic  Statistics  for  Measured  Peak  Flows  Series 

Count 

16 

41 

38 

16 

10 

14 

16 

21 

15 

16 

25 

16 

27 

14 

2 

7 

7 

5 

6 

Mean 

1147 

255 

1037 

792 

1982 

1788 

1503 

3497 

3697 

848 

697 

1207 

170 

640 

1220 

29 

30 

770 

20 

std.Dev.... 

1299 

251 

891 

999 

1337 

1303 

1796 

2493 

3694 

625 

713 

734 

184 

873 

CoefVar... 

1.13 

099 

0.86 

1  26 

0.67 

0.73 

1.19 

0.71 

1.00 

0.74 

1.02 

0.61 

1.08 

1.36 

Skew 

268 

1  88 

1.16 

1.97 

0.77 

0.78 

2.47 

1  80 

2.36 

0.42 

1.73 

023 

1.78 

2.10 

Basic  Statistics  for  Extended  Peak  Flows  Series 


Count 

35 

49 

52 

35 

10 

14 

43 

43 

43 

35 

35 

35 

35 

35 

2 

7 

7 

5 

6 

Mean 

1233 

254 

1240 

834 

1982 

1788 

1414 

3519 

3188 

1115 

850 

1305 

201 

978 

1220 

29 

30 

770 

20 

Std.Dev 

1020 

238 

918 

824 

1337 

1303 

1167 

2050 

2534 

603 

757 

669 

193 

755 

CoefVar 

083 

0  94 

0.74 

0.99 

0.67 

0.73 

0.83 

0.58 

0.79 

0.54 

0.89 

0.51 

0.96 

0.77 

Skew 

2.30 

1.92 

069 

1.61 

0.77 

0.78 

3.31 

1.50 

2.79 

-0  D4 

1.18 

-0.02 

1.28 

0.75 

NFVH  Bui   NFVR  Na/T  NFVR  MM.     NFV*  Spr   Nortfl  Cc. 

Virgin  Vlr 

Gaging  Stations 

La  Varidn   Kanarra    Ash  Cr   So.  Ash  Mil  So.  Ash  fin 

LaadsCr  Santa  Cla/a  Moody  Was  Sovlor  Rv 

Llttta  Cr   Cantor  Cr  Summit  C 

1   Coal  Cr 

Shirts  Cr 

03405420  0940S450  03405490  03405500  03405900  03406000  03406150  03406300  03406500  03406700  03406800  03408000  03408400  03403500  10174500  10241400  10241470  10241600  10242000  10242100 

mmmmTIT]                     2770 

LK 

10600 

K 

i 

5100 

1210 

920 

12000 

502 

1460 

2500 

4360 

770 

4350 

1040 

480 

2610 

982 

390 

5100 

550 

500 

1240 

504 

500 

11000 

1130 

2650 

803 

3400 

1490 

5100 

826 

^■T~ 

3100 

1660 

689 

710 

2770 

643 

1200 

4300 

750 

2600 

536 

39 

00 

4200 

428 

3000 

2370 

3550 

2500 

9000 

30 

00 

2350 

940 

1550 

1980 

1760 

840 

44 

80 

6300 

2910 

25 

40 

1920 

830 

70 

00 

13500 

491 

29 

00 

10000 

835 

1850 

41 

10 

4370 

1500 

580 

1630 

31 

00 

2980 

1400 

1140 

800 

11 

80 

3150 

495 

780 

533 

734 

919 

985 

524 

1250 

13 

70 

1070 

385 

672 

765 

16 

90 

844 

1440 

480 

816 

558 

1700 

290 

366 

208 

35 

20 

2080 

1430 

607 

410 

798 

1400 

488 

408 

11 

00 

1010 

1060 

346 

26 

80 

6620 

317 

211 

833 

2800 

248 

112 

25 

20 

4840 

827 

456 

16 

60 

12900 

455 

1750 

24 

BO 

4690 

410 

1000 

39 

SO 

10600 

532      343 

1570 

956 

2260 

145      254 

214 

668 

1430 

185      626 

651 

22 
*•!                   20 

70 

7410 

396     833 
25     328 

2360 

109 

1 

2573 

4420 

325 

189 

7 

29 

114 

62       8 

1      17 

1210 

i 

2                   4 

is                     58 

172 
6340 

2190 

28 

169 
829 

200 
200 

4 

15 

67 
80 

45      39      259 

■ 
■ 

0      14 
5     27 

477 

144 

13500 

528 

340     1400     816 

2000 

460 

2 

3                   18 

1454 

3100 

107 

210 

130 

5 

3? 

63      546      989 

■ 

1      4 

1040 

370 : 

2 

"j                    43 

4781 

4550 

555 

386 

190 

B 

33 

142      28      285 

I 

2      7 

1600 

273 1 

1 

T]                   13 
jj                   59 

1674 
6764 

4630 
6890 

232 

99 

210 
250 

28 
48 

2980      63      26      628 
2060      81      199      795 

3 

28 

64     24 

1750 

1070! 
144 

I            353      99    2340 

1 

m                                                 16 

1458 

3930 

88 

158 

938 

565      40      791      450 

75      23      22     344 

75 

2 

I                    91 

10057 

22800 

344            1910     890 

2710      776     1810      572 

23      274      70     3340 

440 

1 

r                                             19 

1337 

6840 

800             500      210 

213      49      31      644 

330      300     600     2440 

470 

- 

*                    38 

4434 

13660 

307             872      54 

386      139     1800     975 

12      56      795     4620 

340 

1 
1 

H                    13 

1344 

2660 
2880 

1000             615     580 
600              98 

94      27 
1 290      30 

28 
29 

320 

8      36      122     313 

212 

A                    23 

1803 

236 

152      24     858     278 

790 

2 

7                    22 

2275 

4369 

96             395 

198      115 

481 

1260 

6      25 

8      85 

150 

1 

_J                   19 

1918 

6790 

223             158      510 

84      206 

69 

958 

40      94      355     1820 

220 

17 

40                      7 

192 

3640 

65              52      125 

8      19 

26 

487 

92      109      2 

0     2400 

265 

1 

78                     41 

4195 

11044 

57              16 

80      82 

41 

498 

290      258      37     4440 

526 

1 

74                     12 

1064 

5267 

34             175 

27      48 

IE  ' 

403 

7      36      22     358 

68 

2 

25                      7 

912 

7207 

94             100 

15      28 

28 

487 

6      116      21      70 

486 

95                   13 

1415 

7198 

64             245 

403      132 

52 

528 

9      32      42     570 

501 

2 

04                    16 

2032 

7600 

43             326 

570      175 

62 

1180 

17       19       53      899 

275 

: 

22                   31 
80      128            21 

3664 
2798 

10830 
3650 

100             281 
19              41 

892      139 
43      95 

54 
44 

1210 

16      20      86      515 

640 
431 

406 

3 

67      39     376 

2 

04      242            27 

3454 

9700 

296             113 

278      84 

41 

739 

14 

16      73     651 

358 

2 

72     214            24 

2569 

4740 

606 

802      240 

77 

1380 

8 

57     515     1480 

529 

1 

560            20 
«3                   10 

2091 

4580 

I 

149 
176 

712      36 
12      48 

30 
33 

650 
626 

13 
15 

191      129     695 
28      28     3840 

330 

50      304 

2920 

■                                                     6 

10      346 

1620 

■ 

134 

811       29 

28 

614 

11 

12      20     364 

241 

V                                                 10 
15 

10     1630 
90     950 

7200 
2690 

w 

1 

194 
756 

1170      56 
4420      307 

35 
92 

341 

616 

34 

14 

146      51      840 

440 
344 

74     26 

603 

■                    6 

24      530 

1500 

1  1 

114 

62      20 

26 

181 

6 

i      3 

6      44 

2500 

125 

34      590 

3020 

:/3oH 

116 
138 

125      21 
27      31 

29 

207 

7 

9 

1      15 

1      14 

5       6 
5     35 

140 

237 

74 

M                             ? 

26      756 

1490 

moH 

253 

176 

^            2090     25 

20     2890 

4270 

216 

83      66 

37 

296 

4 

3      5 

4     a 

950 

343 

1 

jj            1429     17 

30     2010 

3920 

398 

3 

58 

147 

56 

1220 

160      188     26 

750 

419 

Basic  Statistics  for  Measured  Peak  Flows  Series 


11 

4 

2 

70 

9 

77 

7 

23 

8 

16 

14 

30 

34 

15 

69 

21 

23 

23 

63 

16 

326 

286 

1760 

2184 

1112 

4325 

1183 

251 

991 

369 

307 

716 

115 

530 

657 

93 

100 

177 

1120 

346 

450 

1619 

833 

3333 

265 

458 

291 

1026 

140 

618 

316 

116 

102 

253 

1047 

257 

1  38 

0  74 

075 

0.80 

1.05 

1.24 

0.95 

1.43 

1.21 

1  16 

0.48 

1  24 

1.02 

1.43 

0.93 

0  74 

3.23 

1.87 

1.20 

1  91 

1  49 

2.61 

1.32 

2.25 

3.40 

1  28 

066 

1  29 

1.30 

1.81 

1.63 

1.67 

Baste  Statistics  for  Extended  Peak  Flows  Series 


35 

4 

2 

84 

35 

84 

7 

34 

8 

35 

14 

35 

35 

39 

69 

35 

35 

35 

63 

35 

230 

286 

1760 

2245 

2365 

5056 

1183 

248 

991 

327 

307 

732 

115 

469 

657 

103 

110 

186 

1120 

359 

262 

1604 

2074 

3846 

219 

348 

291 

962 

138 

410 

316 

102 

91 

218 

1047 

214 

1.14 

0  71 

0.88 

0.76 

0.88 

1.07 

0.95 

1.31 

1.20 

0.87 

0.48 

0.99 

0.83 

1.17 

093 

060 

5.71 

1.67 

1.90 

1.81 

1.79 

2.99 

1.32 

2.28 

3.44 

2.01 

0.66 

1.26 

0.98 

1.76 

1.63 

1.13 

Appendix  III:       Partial  Series  of  Peak  Flows  for  Virgin  River  at  Virgin 


Wat-Yr 

Month 

Day 

Ysar 

Peak  |cfs) 

1910 

1 

1 

10 

2770 

1911 

9 

30 

11 

10600 

1912 

7 

31 

12 

5100 

1913 

10 

27 

12 

12000 

1914 

7 

30 

14 

2500 

1915 

9 

3 

15 

4360 

1916 

7 

26 

16 

4350 

1917 

10 

6 

16 

2610 

1918 

3 

12 

18 

5100 

1919 

9 

3 

19 

1240 

1920 

8 

19 

20 

11000 

1921 

8 

22 

21 

2650 

1922 

8 

31 

22 

3400 

1923 

7 

22 

23 

5100 

1924 

9 

10 

24 

3100 

1925 

8 

25 

25 

1660 

1926 

10 

5 

25 

2770 

1927 

9 

13 

27 

4300 

1928 

10 

31 

27 

2600 

1929 

7 

31 

29 

4200 

1930 

8 

4 

30 

3000 

1931 

11 

17 

30 

3550 

1932 

2 

9 

32 

9000 

1933 

9 

8 

33 

2350 

1934 

7 

28 

34 

1550 

1935 

4 

8 

35 

1760 

1936 

7 

31 

36 

6300 

1937 

5 

8 

37 

1920 

1938 

3 

3 

38 

13500 

1939 

9 

6 

39 

10000 

1940 

9 

17 

40 

4370 

1941 

5 

6 

41 

2980 

1942 

10 

13 

41 

3150 

1943 

3 

9 

43 

919 

1944 

5 

12 

44 

1070 

1945 

5 

3 

45 

844 

1946 

8 

12 

46 

1700 

1947 

10 

29 

46 

2080 

1948 

9 

16 

48 

1400 

1949 

9 

8 

49 

1010 

1950 

7 

8 

50 

6620 

1951 

8 

29 

51 

2800 

1952 

12 

30 

51 

4840 

4 

28 

52 

2230 

5 

4 

52 

2620 

6 

3 

52 

2330 

1953 

8 

1 

53 

12900 

7 

19 

53 

3200 

7 

30 

53 

2080 

1954 

9 

12 

54 

4690 

7 

25 

54 

3280 

8 

4 

54 

1820 

1955 

8 

25 

55 

10600 

10 

8 

54 

4250 

7 

21 

55 

1670 

8 

4 

55 

3220 

8 

8 

55 

2260 

8 

12 

55 

2220 

8 

24 

55 

3390 

1956 

1 

27 

56 

2260 

1957 

6 

10 

57 

1430 

1958 

9 

3 

58 

7410 

11 

3 

57 

3140 

2 

25 

58 

3190 

3 

16 

58 

2190 

4 

23 

58 

1680 

5 

11 

58 

3510 

7 

23 

58 

2780 

9 

12 

58 

5500 

1959 

8 

3 

59 

4420 

8 

18 

59 

2440 

8 

19 

59 

3970 

1960 

9 

1 

60 

2190 

Wat-Yr 

Month 

Day 

Year 

Peak  (eft) 

1961 

9 

17 

61 

13500 

11 

6 

60 

2410 

7 

4 

61 

2800 

8 

3 

61 

2490 

8 

10 

61 

2960 

8 

11 

61 

2490 

8 

22 

61 

2960 

8 

23 

61 

2860 

8 

24 

61 

1920 

9 

8 

61 

2980 

1962 

2 

12 

62 

3100 

2 

9 

62 

2380 

9 

28 

62 

2150 

1963 

9 

18 

63 

4550 

8 

5 

63 

3770 

8 

18 

63 

2110 

8 

19 

63 

1750 

1964 

8 

12 

64 

4630 

10 

18 

63 

2150 

8 

4 

64 

1820 

1965 

9 

5 

65 

6890 

1966 

11 

23 

65 

3930 

11 

25 

65 

1630 

12 

30 

65 

2290 

7 

30 

66 

2190 

1967 

12 

6 

66 

22800 

9 

24 

67 

7220 

1968 

8 

7 

68 

6840 

7 

30 

68 

3210 

8 

2 

68 

1800 

1969 

1 

25 

69 

13660 

1 

21 

69 

2550 

2 

25 

69 

1710 

4 

23 

69 

2110 

5 

2 

69 

1960 

5 

9 

69 

2010 

1970 

8 

18 

70 

2660 

8 

4 

70 

2230 

1971 

8 

21 

71 

2880 

8 

22 

71 

2220 

8 

23 

71 

2070 

1979 

3 

28 

79 

7600 

5 

6 

79 

4260 

5 

20 

79 

3210 

8 

12 

79 

1630 

1980 

9 

10 

80 

10830 

1 

14 

80 

2130 

4 

21 

80 

1730 

1981 

7 

15 

81 

3650 

4 

19 

81 

2110 

1982 

8 

23 

82 

9700 

3 

12 

82 

1990 

1983 

11 

30 

82 

4740 

12 

23 

82 

3290 

3 

14 

83 

2470 

4 

30 

83 

1830 

5 

9 

83 

1820 

5 

28 

83 

2910 

9 

24 

83 

3580 

1984 

9 

10 

84 

4580 

7 

22 

84 

2510 

1985 

7 

19 

85 

2920 

1986 

3 

8 

86 

1620 

1987 

7 

20 

87 

7200 

8 

7 

87 

2550 

1988 

11 

5 

87 

2690 

4 

21 

88 

1810 

8 

1 

88 

1820 

8 

6 

88 

2260 

1989 

7 

28 

89 

1500 

1990 

8 

15 

90 

3020 

9 

23 

90 

1686 

1991 

9 

6 

91 

1490 

1992 

8 

23 

92 

4270 

5 

29 

92 

1790 

1993 

2 

20 

93 

3920 
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Appendix   IV:      Single-site  Frequency  Analysis  of  Annual  Peaks  Series 

Within  the  Study  Area 
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Upper  Valley  Cr  nr  Escalante  -  9336400 
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